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ABSTRACT 


Three analytical models of rotation problems within a 
multi-location, multi-grade personnel system are developed. 
The first model demonstrates the interrelations between 1) 
billet structure, 2) rotation structure and 3) promotion 
structure, and a decomposable transportation problem is used 
to determine an optimal rotation pattern and transfer cost 
while preserving the relations between Dec) ecatlid. Sec heat'o Gs 
In one or more of these three factors are discussed, and the 
model is used to analyze the effects of such changes on the 
total transfer cost. Another model examines the problem of 
transferring personnel between locations and develons a 
parameterization of the personnel selection process when 
transfer costs are taken into consideration. The third model 
examines the problem encountered in Searching for qualified 
people to replace those leaving specified positions. The 
models, although developed independently, are shown to 


interrelate on a macroscopic level. 
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I. INTRODUCTION 


A. BACKGROUND 

In the recent past, concern about the inability of large 
institutions with multi-location, multi-gradée personnel sys- 
tems to transfer personnel and remain within budgetary limi- 
tations has grown rapidly. Transfers within an institution 
are often required to maintain operations of the institution, 
to develop present personnel resources to meet forecasted 
needs of future operations, and to share the burden of 
"hardship" tours in locations of unusual economic, political 
or geographic habitability conditions. On the other hand, 
institutions have become more and more concerned about ef- 
ficient use of their resources. Thus some transfer policies 
have become unodeds but the basic requirements for trans- 
fers remain. 

In order to ostensibly resolve this concern and in an 
effort to develop more efficient transfer policies, several 
rotation models have been developed within the past few 
years. Most of the models so developed can be categorized 
Dye ewomenaractorrsti1es, One ts thessplecifilc intent for use 
of the model; the other is the technique used to perform the 
analysis. Several models have been developed with the intent 
of mechanizing the assignment procedures of such a large 
personnel system, thereby hopefully reducing the workload of 


the personnel managers. Others have attempted to attack the 





more As togss soloed of the development and evaluation of 
alternative transfer policies. The latter may be referred 
to as rotation planning models. 

Several techniques have been used. In the models deal- 
ing with the assignment procedures, network flow algorithms 
have been used by Hatch [17,18] and mathematical programming 
algorithms have been used by Thorpe and Conner [35]. In the 
rotation planning models, several use simulation [27,28,37]. 
Others have used mathematical analytical approaches [8,34] 
and curve-fitting techniques [4]. 

In the development of the analytical planning model in 
chapter II of this thesis, the interdependence of the per- 
sonnel inventory SCT Weuu7 cemcalleds them bil let#S tnueture 
PiImoucwoute uve aire sh che the promoticn Structure ana once rova- 
tion structure of the institution become quite apparent. 
Several models of personnel inventory structures have been 
developed. Among these, Markovian schemes have been used 
(1,29,30] as well as entity-simulation schemes [19,33]. More 
recently, Hayne [21] has developed an approach to Markovian 
schemes which allows two characteristics to be used to de- 
scribe the states of the Markov model. Also, Marshall [26] 
has compared the use of a Markov model with the use of a 
cohort model for use in the prediction of the future person- 
nel inventory of an institution. Grinold and Marshall [16] 
provide a thorough development of cross-sectional and longi- 
tudinal models, along with discussions about optimization 


and data formatting. 
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Only one model has been found that explicitly attempts 
to interrelate promotions, withdrawals, reracions wand inven- 
tory. The model was developed by Hatch [20] and uses the 
entity-simulation technique with little documented analysis. 
The model in chapter II gives an analytical approach to the 
problem, and further develops an optimization scheme to 
evaluate different transfer policies. 

The model presented in chapter [II is concerned with the 
assignment process, so may be classed among those models. 
However, the emphasis is not on mechanizing the assignment 
process but rather on presenting different policies that the 
personnel managers might use to guide them in the process. 
PI rOreeiseniadesno relate the efrect of Che policies on the 
expenditures required to maintain the transfer process. These 
models provide planning aids which take into account factors 
previously ignored in the planning precess. The model of 
chapter II is designed to be used on an interactive basis 
by the planner. 

Lieis noted that “the majority sor the reterences listed 
above are to unpublished technical reports which were in- 
tended to focus on accounting problems and not on in-depth 
analysis of implications of policy changes. Although this 
thesis deals with the same or similar problems to those 
found in these reports, the emphasis here is on analytical 
approaches which give insight into the complex interactions 
among retention, promotion, rotation and inventory. 

Other Prey tical work has been done in manpower planning 


and forecasting. Among early workers were Bartholomew [2], 


1] 





— 


Blumen, Kogen and McCarthy [3], Young and Almond [40] and 
Gani [11]. This work has focused on forecasting of inven- 
tory and labor turnover problems. More recently, Charnes 
and Cooper together with other authors including Niehaus, 
Sholtz, Stedry and Klingman (see list of references to [6]) 
have worked extensively in applications of mathematical 
programming to the mechanization of the personnel assignment 
process with emphasis on civilian manpower in the Department 


of the Navy. 


B. CONTENTS AND SUMMARY 

In this thesis, three analytical models which describe 
rotations in multi-location, multi-grade personnel systems 
are developed. Such systems are the basis of the personnel 
Structure of many large institutions. The U.S. Navy is an 
example. 

The three models are essentially independent, i.e. none 
depend on either of the other models for input or output. 
However, they are shown to interrelate macroscopically. 

That is, the insight provided by one may change the manner 
in which another is used. 

In chapter II, some basic balance equations which relate 
the promotion structure, the rotation structure and the bil- 
Retestructure of a llarce institutioneare discussed. Two 
basic equations are developed in section B. These relations 
tenmusecd repeatedly throughout the remainder of the chapter. 
Beectilou se COntains the development of three different ways 


in which the billet structure may be characterized. These 
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prove useful later in the chapter. Sections D, E and F con- 
tain the development of the methodology required to demon- 
Strate the effect of changes in the personnel system on the 
basic balance equations developed in section B. In section 
D, the necessary calculations to show the effect of various 
changes in the billet structure on the promotion structure 
holding other factors constant are performed. In section E, 
Ener et tcet of changes Wm the rotation structure on the ypro- 
motion structure, holding the billet structure fixed, is 
described. Finally, in section F, changes which maintain a 
fixed promotion structure are described. 

Thesiext threceysections of chapter 1] aré=dévoted to the 
development of a methodology to compare the relative effec- 
tiveness of various changes in the system. Basically, from 
sections D, E and F we find that the same change in one 
component of the system can be compensated for by various 
changes in the other parts of the sysitem. The choice of an 
"optimal'' way of compensating for a specified change is then 
of interest. In section G, the results of the chapter up to 
that point are summarized, then the problem of comparative 
analyses of alternatives is introduced. A transportation 
problem format is mecdueed in section H as a method of 
errecting the comparision. Sensitivity tanalysis on the 
Pransporedtion problem is discussed¥in section 1, both in 
meneral and for a specific case - that of a change in the 
miebet Structure for which the promotion structure is to be 


compensatory. 
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In section J some of the assumptions that were used 
earlier in the chapter are relaxed and the limitations that 
result with the relaxations are demonstrated. The chapter 
concludes with a discussion of possible extensions for the 
model. 

In chapter III, attention is shifted from the flows a 
personnel within the institution to the problem of selecting 
a suitable replacement for an individual who is leaving some 
position within the institution. The model that is developed 
is called an Ehrenfest Decision Model because of some basic 
Similarities to the Ehrenfest Urn model. 

Following an introduction, in section B the Ehrenfest 
Decision Model is developed in the special case where the 
institution is confined to two separate locations. The 
model is developed under a simple decision rule, and then 
the decision rule is generalized. Approximations prove use- 
ful and are developed. In section C, the model is extended 
to institutions with three locations. Some comments are 
made. about the difficulties to be expected if more than 
Pitecem10Catlonssare to be described, 

The last section in the body of chapter III contains the 
development of diffusion approximations to the model of sec- 
tion B.| Results from this approximation are useful in the 
description of transient cases. The conclusion to chapter 
IilI contains some remarks on the macroscopic influence of 
the results of chapter II on the interpretation of the re- 


sults of the Ehrenfest Decision Model. 
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Our attention in chapter IV is on the problems encoun- 
tered by the manager who is charged with making the selec- 
tion of a suitable replacement. It is seen that a "good" 
choice of model parameters from chapter III may not be as 
good as was thought. The chapter begins by assuming in 
section B that the manager is operating with no restrictions 
on the number of personnel from which the selection is to 
be made and that he searches until the "right" replacement 
is found. The description is then extended to the case when 
the manager has a limited effort to be applied to finding a 
replacement, discussed in section C. In section D, restric- 
tions on the number of personnel from which the replacement 


may be found are introduced. Restrictions are then made on 
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mach the numcos 
being the basis for section E. In the conclusions to the 
chapter some possible extensions and also the influence of 
the results of chapter II on the results of this chapter are 


discussed. 
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II. RELATIONS BETWEEN PROMOTIONS, ROTATION RATES AND BILLET 
STRUCTURE IN A CLOSED INSTITUTION 
A. GENERAL DISCUSSION 

In this chapter we consider a manpower system for an 
institution whose personnel are in various locations and are 
further classified into ordered grades. An example of such 
an institutional manpower system is the enlisted force of 
the U.S. Navy where the locations could be the various ports 
and shore stations and the grades could be the nine pay 
grades within each skill category, or rating. 

Personnel are rotated among the various locations; they 
are promoted between the various pay grades; requirements 
for a given skiil category and pay grade can vary with both 
location and time period. We develop relations between pro- 
motions, rotation rates and the billet structure of such an 
institution. These relations are used to analyze the ef- 
fects of changes in these factors under various assumptions. 

Consider an institution in which personnel are classified 
into hierarchically ordered grades and in which personnel are 
hired only into the lowest grade. Again the enlisted force 
of a military institution serves as an example. 

In this text, An;k denotes the kth assumption in chapter 


m; Dn;k denotes the kth definition in chapter n. 


DZ2;1: A personnel system which hires only into the lowest 


grade is called a closed system. 
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The legacy of past recruiting, training and advancement 
severely constrains the range of current feasible decisions 
which can be used to control tke personnel structure of an 
institution with a closed system. In addition to maintain- 
ing operation of the institution, the personnel managers 
must also ensure the development of present personnel re- 
sources to meet forecasted needs of future operations. Thus 
they assume responsibility for the training and advancement 
sequence to meet future requirements. 

An institution using a closed system may wish to rotate 
personnel periodically. There may be several reasons for 
the rotations. For instance, one purpose for periodic ro- 
tations might be the desire to have the personnel exposed 
to a variety of job types. An ability to perform well in a 
variety of positions might be a measure used in the insti- 
tution's promotion scheme. In this case, rotations might 
prove to be "career enhancing" for the individual. 

For geographic, political, economic or other reasons, 
some of the locations where the institution uses personnel 
may be less desirable than other locations. The institution 
may desire to establish an equitable rotation plan influenced 
By such location factors. Specific positions may require 
extensive periods of family separation and this may be rec- 
ognized by a policy of follow-on tours in positions more 
conducive to the maintenance of a close family life. 

PnhePinstitution may find that the rotation policy af- 


Be€tsceitS retention. In a closed system the retention of 


UY ) 





highly capable personnel is a significant factor in maintain- 
ing the caliber of the future leaders of the institution. 

For instance, institutions involved in the extraction of 
raw mineral resources may conduct operations in areas deemed 
inadequate for continued habitation. Government employees 
are often utilized at "undesirable" locations. Without an 
effective rotation plan, both the mining and governmental 
institutions may have poor retention and hence possibly less 
effective future operations than would be possible with a 
more effective plan. . 

Whatever the reasons for establishing the rotation 
policy, the institution might choose to set minimum and 
maximum tour lengths for various positions at the various 


Weoesticns. Minimum tour lengths might be sct to ensure thac 


r 


an individual's performance in a specific position can be 
properly evaluated. Maximum tour lengths might be prescribed 
in response to any combination of the above hypothesized 
reasons. Other factors peculiar to a particular institution 
might also be considered. In this thesis we assume that ro- 
tation planning is based on a set of given tour lengths. 

Let the tour length associated with a position at location 1 


with grade k be denoted T(1;k). 


tae ee the rotationsratestor a location"! andga grade k, 
fenceed Rii;,k), is the reciprocal of the tour Iength for a 
position of grade k at location i, i.e., R(i;k) = eelae 
UC Sivebetestructure ot Air Institution 1S simply the 


number of positions of each grade at each location. 
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D2;3: The biliets for a grade k at a location i, denoted 
b(i;k), give the number of personnel required in grade k 
at location 1. 

The personnel rotation scheme to be investigated is 
essentially based on close examination of the billet struc- 
minc. mens meen by 2Z;35, Ehere are two attributes which each 
individual in the system is assumed to possess. The first 
attribute is location. The second attribute might be the 
pay grade of the individual, the skill group to which the 
individual belongs or another characteristic which is of 
interest to the rotation planners. 

In the next section, we develop two basic relations 
which are used repeatedly throughout the remainder of the 
chapter. One arises from basic flow considerations, while 
the other comes from the assumption that the institution 
being investigated has a closed personnel system. The rela- 
meronssretlect the "Interactions between the billet structure 
uecile S/S tCiiymeme rotation Structure Gf the system and thie 
promotion scheme used in the system. 

Section C contains the development of three different 
ways in which the billet structure may be charactcrized. 
These prove aoc Wiel anen san the thesis. 

rimscCELOlS@D ss andere wesdevelop the methodology re- 
quired to demonstrate the effect of changes in the system 
pre ercebaste relations developed in ‘section B. The next 
meiroce sections, G, H and 1, are devcted to the development 
of a methodology to compare the effectiveness of various 


Changes to the system. 
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In section J, the effect of the relaxation of some of 
the assumptions in section B on the tractability of the 
problem is discussed, We then conclude the chapter with a 


discussion of some possible extensions to the work presented. 


B. DEVELOPMENT OF THE BASIC BALANCE EQUATIONS 

In this section, we develop two basic relations which 
are used repeatedly throughout the remainder of the chapter. 
One arises from basic flow considerations, while the other 
comes from the assumption that the institution being inves- 
tigated has a closed personnel system. 

It is assumed that the rotation rates R(i;k) are given. 
It is also assumed that the billet structure is known. By 
scanning the current list of employees and comparing each 
individual's date of assignment to his current position with 
the stated tour length, the rotational needs within the plan- 
ning horizon can be determined, assuming that every vacated 
position must be filled. For Simplicity the planning horizon 


is now taken to be one year. 


eee ne TOtavwonas nieedSewi thineplanning period n, de- 
iocedema( i kon jesane the requirements ate location 1 for per- 


sonnel who have second attribute k. 


For ease of discussion in the following, the second at- 
Peepeieee ls taken to Despay wyade. | AS Stated ingthe previous 
section, ‘the attribute used might also be skill group or 


elemmr daetan of antéeérest) to the plannegs: 
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D2;5: The number of personnel of grade k at location i 
whose tour is completed during period n (and hence must be 
rotated to another location) is denoted a(i;k3;n) and called 


the availabilities at location i of grade k during period n. 


Note that the locations are not necessarily physically 
separated. For instance, the U.S. Navy may designate one 
location as San Diego-Sea duty and another location San 
Diego-Shore duty since its rotational needs require some 
distinction between the two types of duty. 

Since availabilities arise through completion of tours, 
there exists some transformation on past requirements which 


provides the current availabilities. 


A2;1: There exists a square matrix Q which provides the 
transformation on past requirements giving current avail- 
abilities. Further, denoting by q(m,k) the element in the 


mth row and kth column of Q, we assume 


a(i;k3n) = in r(i;mj;n-T(i;m))q(m,k), (2.1) 


mor abl locations 1, grades k and periods m of interest. 
Note that q(m,k) gives the fraction of personnel who 
sono letomthelr etour ein grade k .meiven theysemter a location 
my grade m. Thus if there are G grades, Q is a G-by-G 
aver 1 X | 
The existence of Q postulated by A2;1 implies that a 
linear transfcrmation on past requirements is sufficient to 


determine current availabilities. The possible lack of 
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agreement of this assumption with the structure of the sys- 
tem under investigation must be weighed against the tracta- 
bility of the proposed models The purpose of the curnent 
investigation is to aid in the development of insight into 
Lake system. As such, it is seen as a planning tool rather 
than a rigorous forecasting model. In this light, the trac- 
tability and resultant analytical results seem quite helpful. 
It should be noted that the matrix Q has the structure 
of a finite state transient Markov chain. The significant 
difference between the Markovian model presented herein and 
Markov models used elsewhere lies in the interpretation of 
the matrix Q. In most cases, for example Bartholemew [2], 
Blumen SR a [3], Charnes, Cooper and Niehaus [5], and Young 


139}, tleé Markey transition matrix is defined to provide the 
probability of transition between states during the course 
of a fixed time period. The period varies, sometimes being 
one’ year, sometimes one or more quarters of the year, etc, 
but the transition matrix is assumed to describe the flow 
murine a period. 

In this thesis, we assume the transition matrix to be 
defined over tours, and have established that tours vary over 
grades and locations. The resulting matrix of transition 
fractions (since we have defined fractional flows and not 
probabilisitic transhtions)) 1s thus)different from those of 
previous Markov models since there 1s no fixed period for 
which the transitions are measured. 

Thus the elements of the matrix Q give the transformation 


on grades confounded by possible variances among tours. If 
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all tours within the institution are the same length, the 
matrix is interpreted in the same manner as previous models. 
Suppose also that there exists a 1xG vector, w, where 
w(k) gives the fraction of the requirements for personnel of 
grade k which withdraw from the system prior to or during 


Enempcriod) Of their meéxt availability. 
A2;2: The withdrawal vector, w, is strictly positive. 


DUC A cae wi race homo 7 equi rementcmmonrs personne I 
of grade k which will be available upon completion of a tour 
will always be less than one. It is noted that A2;2 is more 
rigorous than required mathematically for the succeeding 
development. However, observation of the Navy personnel sys- 
tem indicates that AZ;2 holds for that system. If a partic- 
ular institution does not exhibit the withdrawal pattern of 
A23;2, then for all grades g from which there is no withdrawal 
one may set w(g) = € where € is a small positive real number. 

ie each ellémentimwoswO- Gg (lM. (ia) useimterpreted=to give the 
meactlOnmer penSennelowmo Entered their present tour (in 
response to a requirement) as grade m who will be available 
for reassignment as grade k upon completion of the tour, 
miem the .sum over#k of q{(m,k) will be lics& than one. In 
act, a q(m,k) + w (m) ol t:Orecel Grades. me— lt... ,G. 

If the institution is neither gaining in population nor 
decreasing in population, cvery withdrawal must be matched 
Mya FeCruit. Simce the systém is assumed to be clos®€d, all 


Meenults must enter in grade 1, the lowest grade. Then the 
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number of recruits during period n, denoted r,@, 1S given 


by 
LC 
ro(m) = 52) 42, riskin-TG35k))w(k). (2.2) 


Equation (2.2) simply relates the recruiting requirements 
to the assumed withdrawal rate over all locations i, 
ieee. Geamad grades k, k= 1,...,G and to historical 
requirements. 

Let us assume that the requirements at n-T(i;k) are 
Pdutemrormtne requirements in perlod nm for all Locations 1 
and grades k. This is a “weak'' form of stationarity in that 
some cyclical variations may occur so long as the require- 
Mentsmhepcau cvery I (iskjaperiods. 

Let a(ij;n) be a 1xG vector giving the availabilities at 
focatlonei during period n, where the e@ements, of afi;n) 
cman ik od ilben ive |Ctmm (a. nje besa IxG veckor Sivine 
Liem cgulrementseat location 2. Lhem equation (2. ljm@may be 


written for any period n as 


a eter ( 1 sh Oe (274 9) 


If all those available to move must move due to the con- 
Siderations of the preceeding section and if all requirements 
of the system must be met, then 

L 

i2y a(ijn) = 521 1) 31) ae (r,(n),0,...,0), 
eiere the last vector on thétright-hand-side gives the in- 
Crease in availabilities of personnel of grade 1 duc to 


recruitments. 
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SUDStILULIneg elation (2Z.1') into the left-hand-side 


of the above equation and rearranging terms, we obtain 


L 
(I-Q) i241 r(15n) a (r,(m),0,...,0), 


1 


or, if (I-Q) ~ exists where I is the identity matrix of 


appropriate dimension, 


1 Velie ot) )) = r,(@)N,, (2.5) 


where Ny is the first row of iO. The notation N = ene 
follows Kemeny and Snell, [23], who call N the fundamental 
matrix. Since Q may be interpreted as a substochastic matrix 
due to A2;2, the fact that N exists is well known. Kemeny 
and Snell ([23], p. 46, Theorem 3.2.4) give the interpreta- 
MoM Oluethce (init clementeot Neas the €xpecté@ number or 
Pioluseto state 7 Dy an drbitrary entry into sitate®i. Each 
element of the vector equation (2.3) may be stated -Eyr(isksn) 
= r,(a)N, (Kk). Thus, element by element (2.3) equates the 
total requirements for personnel of grade k during period n 
MemenCeproduct of the number of entries inte the*®system dur- 
mre period mn and the expected number of vilsiits)¥to®grade k by 
an arbitrary entry into the system. 

- A disadvantage of the above model formulation is that 
ene factors under the control of management personnel such as 
PielLGcty structure and rotation raites are not explicit in the 


eeructure of the relations (2.2) and (2.3). Thus, the ef- 


fect of varying one’ or both of these factors is not obvious. 
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We now introduce the planning factors explicitly into 
relations (2.2) and (2.3). Due to the assumed cyclical na- 
fure or the requirements, “1.e.—~ens assumption that r(i;k;n) 
= r(i;k3;n-T(i;k)) for all locations i and -grades k, every 
Sic Cie thesp (13K) billets for wemsionne!l of gradewk at loca- 
tion i must be replaced every T(i;k) periods, and during 


this time each will be replaced only once. Hence 


cl) gue 


b(isk) = 7g 


pi fen), 


where the period in which the counting starts, m, is arbi- 
Eratyes it thessystemyis stationary inythe usual sense, mee., 
r(i;k;n) is independent of the period n, the above summation 


gives b(i3;k) = r(€i;k3;n)T(i;k). 


A2;3: The requirements r(i;k3;n) are independent of the 


period n. 


We now suppress the n in the previous notation, setting 
mit;ki:n) = r(isk) for all locations i and grades k. Now 
requirements r(i;k) are given by b(i;k)/T(i;k) or since 
R(izk) = T(i3k) +, r(i:k) = b(i;k)R(i;k). Substituting into 
(2.2) and (2.3) respectively, we have 


b(isk)R(i:k)w(k) (2.4) 


L 6G 
ope” ea Ate 


and 
Neawie): ; 
Joh UE ae 


1 Diipsigii( cielo «ee lene, «4G. 
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The last equation may be stated more succinctly as 


rN, = 1B OR (ZS) 


where the symbol © means element by element multiplication 
ot Lie smatricesmBeand) R and i is! the stumjvector of the jap- 
propriate length. Equations (2.4) and (2.5) are called the 
"basic balance equations" for the system. 

Much notation has been introduced in this section that 
is used extensively in the sections which follow. For this 
reason, a summary of notation is included in Table 2.1. 

The next section contains a discussion of some alterna- 
miVve woys inewhich the buliet #stnucture may be characterized. 
These are of use in subsequent examination of the implica- 


tions of the basic balance equations, (2.4) and (2.5). 


C. FRACTIONAL BILLET STRUCTURES 

.This section develops three different ways in which the 
Daublet siructure*may be characterized. The characterizations 
are not independent, but each is useful in specific situa- 
PlOousmiLater in the thesis. 

her epe DCm LMCmLOtd! GUMben Of pemSonnc lum Liecesystem 


Ie G i 
1.e. b = i234 Ke bm) > Then the?(1,k)th element of b 4 


B 
gives the fraction of personnel whicheare of Grade: keat /16- 
Pc LOwm@etem HenOtemtinis fractional distribution matrix by IT, 
with (1,k)th element y(i;k) = b(i;k)/b. One sees immedi- 

ately that the fractionalematrix can be substituted for the 


PeelctumStTucture matrix B in (2.44) and (2.5), and that these 


Paiattons dre linear in the total system size b. Hence 
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Tab leeZ. i. 
SUMMARY OF NOTATION 


L is the number of locations within the system, indexed i. 


Gis the number of grades (or other second attribute) 
Within the system, indexed k. 


T(i3;k) denotes the Tour length in integral periods for 
a person of grade k at location i. 


R(i;k) denotes the rotation rate for a position of grade 
k at location i and equals 1/T(1i;k). 


Ris a LxG matrix with elements R(i;k). 


b(1;k) denotes the number of billets of grade k at 
location 1. 


B is a LxG matrix with elements b(i;k). 


r(i;k3;n) denotes the number of personnel of grade k 
required to be assigned to location i during period n. 
In the stationary system, the n is suppressed. 


Terror amas a LxGevectonm with welementsS Fr(i;ken) or 
rir, kee respectively. 


hiMesOtor,sd.ixGhevector defined by re= (r(1),...,r(L))-. 


a(i;k;n) denotes the number of personnel of grade k 
which. are available to move from location i during period 
een thesstattomary sgsitem, the m 1s ysuppressed.~ This 
1s an element of a(i3;n), a(i}, a(n) and a in the same 
manner as r(i;k;n) is an element of the corresponding 
requirement vectors. 


Q is a GxG matrix giving the fractional flows among 
grades during a tour. 


wis a 1xG vector giving the fraction of personnel who 
enter a tour in grade k who choose to leave the system 
upon completion of the tour. Elements referred to as 
withdrawal rates. 


ee COG r,@)), the number of ré€cruits (during period n). 


O 
Ny; telwomrcCeros whach 1s the first row of ety > 
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policy alternatives may be examined using the fractional 
distribution matrix and results applied to an institution 
of any size. It is noted that [ has dimension LxG. 

Let b(+;k) = 7 b(i;k), the total number of personnel 
of grade k in the system. Denote by © a GxL matrix with 
fie 1) th elementwd(k:i)m=-8b (i 3k) /b(t+ 3k). That is, the 
ike) th element of © pives the fraction ofjpersonnel of 


grade k who are located at i. Note that 


y(isk) = $(k;1)b(+;k)/b. Cre 


SinWariv. let "b{1s+) = b(i;k). Thus b(i;+) gives 


G 
k21 
the population at location i. Denote by 9 a LxG matrix 
With elements 9(i;k) = b(i;k)/b(i;+). Here the (1i1,k)th 


element of 9 gives the fraction of personnel at location i 


who are of grade k. Note here that 
y(isk) = 8(i;k)b(i,+)/b. (Ze 7) 
From (2.6) and (2.7), we see immediately that 


b(i;k) 


b(+3;k) 6(k3i1) (2. Sa) 


bi s+ 0 (1 : ke (2.8b) 


and that the row sums of 6 and 9 are unity. Summing (2.8a) 


over k and (2.8b) over i gives 


b(i3;t) b(+3k) $(k31) 02,92!) 


a G 
2 


k=] 


and 


b(+5;k) Dilan) 0 (1 ke (2.9b) 


L 
Bs) 
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Define Dh: as a 1xL vector with elements b(i;+) and 
beg as a 1xG vector with elements b(+;k). Then the pre- 
ceeding relationships may be expressed in vector-matrix 
notation as bys = b.é and b.¢ = b, +9. Substituting one 


into the other gives us the relations 


b 


L° b, 0 (2enE0 a ) 


and 


b.p = b., 90. (2.10b) 


Lemma 2.1. The matrices 0% and $0 are non-negative and have 
row sums equal unity. 
Proof: Non-negativity follows from the definition of 0 and 


®. First consider 90%, Then 


I 


mz1 °C) »™) 52, 9,1) 


G | 
mzy SC) »m) = 1. 


The result for #0 follows by a symmetrical argument. 

By virtue of the above lemma, both 06 and 60 have the 
properties of a Markov matrix, which we assume will be 
regular. Therefore, using the well known steady state dis- 
tributional result of Markov theory (for example, see Kemeny 
and Snell [23]) equations (2.10a) and (2.10b) have solutions 
by. and big respectively which mire unique up to’ a_multipli- 


Sacive. constant. Hence given the. matrices ® and 0 and the 
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total system size b, the marginal billet structures by loca- 
tion and grade are uniqucly determined. 

In summary we have characterized the billet structure 
in three forms, [, ©, and 0. These are used in subsequent 
sections in the analysis of system changes brought about by 


changes in the billet structure or rotation rates. 


D. ABSORBING BILLET CHANGES THROUGH CHANGES IN THE MATRIX Q 
1. Changes in a Single Element of B 

Suppose the population of personnel of a specific 
grade at a specific location is increased or decreased. 
Suppose also that there are no other changes in billets 
mirouchout the system. Thus the total population of the 
institution changes by the same amount. In order to main- 
tain the basic balance equations ((2.4) and (2.5)) there 
have to be changes in some other factors. In this section 
we fix tour lengths and withdrawal rates and observe the 
effect of the billet change on the matrix Q. If we assume 
the second attribute which describes the population is pay 
grade, then Q is a description of the promotion scheme. 

Suppose we wish to change the number of positions 
ae location & in grade g, i.e. the quantity b(4;¢8). “We 
wish to express ry and thesvicic lor Ny as fumetions of the 
ehange in b(%;g). To do this, we parameterize the popula- 
tion at location 2% of grade g by denoting the population 
following a change as b'(2;g) and setting b'(2£;¢9)=(ita)b(2;g) 
eee eOe1S tne fractional change in b(2;¢e), the control var- 


table. It ws straightforward to show that ry is linear in 
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a. From (2.4) we obtain 


ro(a) = ry + w(g)R(2sg)b(23g)a (2.11) 
where 
iF G 
a 324 Kz] BCI; kKIR(: kK) whee 
Alc Oee tron 2.aeand (2.11), 2t 15 trivial to show 
re p (k) 
N, (k30) r+ Ww AWC eS. * k#g (2.12a) 

where 


L 
p(k) = ,2, bG3kK)RG;K), 


or p(k) = r ON, CK) from equation (2.5). We also obtain when 


¢ 


k=g 


-y) = (gs) + RO g)b(23g)a 
Nj (30) = 2S RUase)b UB WCB) oie 


Recall that N 1 


is the first row of (I-Q). Thus 


1 
the above equations do not show explicitly how Q changes 
with a. In order to proceed we assume the following special 


peructure of FQ, 


A2;4: The promotion scheme used by the institution allows 
no demotions; no individual is allowed to advance more than 


one grade during any given tour. 


Reca Piethetethe elementeq (mj) of QO is the fraction 
of personnel who complete their tour in grade k, given 


mieyeentercd their present tour in. grade m. For the 
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special structure given by A2;4 let q(m,m) = q(m) and 
q(m,m+1) = p(m). All other q(m,k) are taken to be zero. 


Hence Q is seen to haver*tne following structure: 


p(1) 0 0 0 0 
q(2) p(z2) 0 O 0 
0 q(3) p(3) 0 0 


0 0 0 q({G-1) p(G-1) 





0 oO 0 0 q (Gees 


In this special case, let n{(i) = (lecqin a Then Ny can 


be shown to have the following structure: 


Ge 
Ny = (n(1), n()p(1)n(2),---.n(G) pl, n(k)p(k)). (2-13) 


Nien chismrorm and recalling) that q(1) + p({i) + w{i)"=81, we 
can iteratively calculate the elements of Q from Ny for any 
fixed w. We note that with the special structure imposed 
by assumption A2;4, the Irecuercation xe ele (al ews Saye 
murce 5S, the gwotation rates R and the withdrawal rates WwW, US 
Sutficient to determine a unique promotion matrix Q. 

We now give an example, the basic elements of which 


are used throughout the nest of the chapter. 


Example 2.1. Suppose we have an institution with four loca- 
Proms andg§tive pay grades. Suppose the billet structure B, 
the rotation rates R and the withdrawal rates w are given 
by: 
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300 240 
600 455 
a 300 240 
1140 600 
BRS Jes 
Bae SBR 
sy 1 1 
Foes 5 


omidmwe=@ (Wels 1.2 55 24). 


180 70 
D500 Pet 0 
180 120 
440 280 
e353 _ 
wid ae, 

1 


35 


60 | : 
80 | 


Then the promotion matrix Q required to maintain 


the given billets and tour pattern is 


~354 .546 
[ . 306 
0 = 
Example 2.2. Suppose now 


at location 4 in grade 1, 


394 
BAT Lo 529 


etOU oe 
7G 


iIgil. 
00° 


that we wish to change the billets 


J. Cam bi 4; 


, to 800 from the cur- 


rent 1140. Maintaining R and w fixed, we ask for the changes 


in Q required by this change in billet structure. The new 


O1Ss 


a O72 eos 
. 306 


. 594 
ayy 2 329 
.489 
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Nope ent thesdecrease in billets of grade 1 Has 
resulted in an increase in promotion rate from grade 1 to 
grade 2. This is to be expected, since a smaller community, 
i.e. the new population of grade 1, is being required to 
support the same population of grade 2 as was required pre- 
viously. Hence a larger proportion of grade 1 will be pro- 
moted during each tour. Perhaps the surprising result from 
example 2 is that only the elements q(1) and p(1) changed 
with the®change in b(4;1). In fact, it is true in general 
that for a change in billets of grade k, only those elements 
of Q, q(i) and p(i), for i less than or equal to k, will 
change. We now show this result analytically, but first 


we need the following lemma. 


Lemma 2.2. Assumption A2;4 and equation (2.13) imply that 


the element q(k) of Q may be expressed 


k-1 
q(k) = 1 - NEY (1 - 52) w(3)N,C) }. (2.14) 


Proof. From (2.13) we have 


I 


N,(1) = n(1) = 1/(1-q(1)), 


from which 


1 


q(1) = 1 - Ny (1)" 
Assume now that (2.14) holds for k. We show that this im- 
Mmonese tivato(2.14) holds for k+l. First we observe that 
(2.13) implies that 


Ny CK) PK) 
eee Toagtet) 
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From the above, we obtain 


Ni (k) p(k) 
q (k+1) =e Waly 


But we also have that p(k) = 1 - q(k) - w(k). Substituting 
for p(k) in the preceeding and cancelling terms, we obtain 


the result. Hence the lemma is shown by induction. 


Theorem 2.3. For a fractional change a in billets of grade 


g at location &, the elements of Q, under A2;4, are given 


by 
TK {C,(k) - w(g)R(2sg)b(2sg)al, k=1,...,g-1 
ry > Ey w5)e(5) 
q(k3a) = ie WS) 0 (8) Fi R(2sg)b(23e)a ’ k=g 
C, (k) me kot te G- 
where 


k-] 
C,(k) = p(k) - x, + 21 W(j)e(j). 


Broot. Ihe proot tollows by straightfonward substitution 


of equations (2.12) into (2.14). 


Example 2.3. Suppose now that we wish te change the billet 
Pemucturec. otmexample 2 by Setting b(1;5), the billets in 
Hoecation 1 of grade 5, to 1000. Again fixing R and w, we 
Hck tor the promotion matrix required to maintain the spc- 


cified structure. We obtain 
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Lode oo 
056 .644 
Qs »156 .664 
UWA, 5 [ ae4 
.600 ° 

Example 2.3 illustrates that there are limitations 
on the ability of the system to absorb changes in billets 
only through changes in the promotion scheme used. The 
matrix Q must be non-negative; hence the element q(4,4) 
indicates that the promotion scheme will not support the 
desired change in billet structure. 

Next we develop relations which provide for a change 
in a specific billet to be absorbed throughout the institu- 
tion in such a manner as to preserve both the total popula- 


tion and the population within each grade. 


Z. Change in One Element of B Absorbed Proportionately 
Throughout the System 


Suppose now that the increase or decrease in the 
population b(23;g) is restricted so that neither the popula- 
tion in grade g, b(+3;g) nor the system population b is al- 
lowed to change. Thus the location populations D(i;+) must 
change. There are two cases that would be of general in- 
merest in this instance. 

In the first case, a sub-unit of location & may move 
aoe aesunit to location k, k#2. This would be the case, for 
instance, if a firm engaged in heavy construction completed 
a project in Idaho and simultancously acquired a contract 
mie et zoids) Ihis case could be easily handled by two itera- 
tions of the scheme in the above sub-section. 
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In the second case, the firm might complete opera- 
tions in one area and not immediately desire to move every- 
Bneeto a especific new area Of operations. In this case, in 
order to retain valuable expertise the firm may redistribute 
the excess personnel throughout the system. We examine this 
case, assuming the distribution is to be made so that each 
location i gains in population in ratio to its population 
of grade g. Now let b'(23;g) denote the population D(%3g) 
following the change, where b'(23g) = (lta)b(23g) = 
pity) (lta) ato;2)o9 We know that the sum over 1 of 79(931) 
is one. Also, letting ¢'(g;2) = (lta)¢o(g3;2) and 6'(g3;1) = 
(1+8)o(g3;1) for i#% to ensure proportional distribution, we 
must ieReoommet 3 the sum over i of 6" (234) is one. We note 
that the Dare notation implies the value foliiowing the 
change. The accumulated changes must sum to zero. Hence 


B i52 Oe Jee agg sf) = 0 


Recognizing that ide o(g31) = 1 - o(g3;2), we obtain 
B = - ad(g5;2)/(1-o(g;2%)). 


We now have that the new billet size for personnel of grade 
feat 1OCatlons 1, 27°, 1S)given by "b'{1;¢) = b(+;3¢) (1+8) o(g 51) 
-= (1+8)b(i3g). We are now able to develop the expressions 


ferethe total wecruits ie and Cie mC ClO man and thus Q, from 


dis 


the basic balance equations. By straightforward substitution 


and some algebra, we obtain 


r Co) eS Kiw(g)a, (21'S) 


Bye 





where 
K, = b(Qsg)R(t3g) - pe 
1 a8 8 1 - eee: 


We also obtain, for k#g, 


Die, oc) Rh), 


“iM c= 


1=1 
172 


p (k) 


Ny (k 3a) = r+ K.w(g)a (2.16a) 
O Ht 
and 
Pomp (oe Kye 
N, (g30) r+ Kywg)a ’ (2 ob) 


where 9 is given in (2.12). 

We note that (2.15) and (2.16) correspond very closely 
to equations (2.11) and (2.12). In fact, the only difference 
is that in the relations developed in this section, the con- 
stant Ky takes the place of the constant product in (LAER en) Rare 
of the previous work. Hence, the results of theorem 2.3 ex- 
tend to this case with the proper change of constants. 


We also note that Ky may be written 


Sy Saar (b(+5;g)R(23g) - vr CO)N, (g30) ). 


The interpretation of this equation must await future dis- 
cussion, but we can note that should the term in parentheses 
be zero, then equations (2.15) and (2.16) would be constant 
Tl oO. 

With assumption A2;4 we can now obtain the new pro- 
motion matrix O using (2.13) and (2.14). *We conclude this 


aibescctrom with anegexample. In the? neéxt?sub-section we 


oe, 





extend our results to the case where the entire force at 


some location is scaled up or down. 


Example 2.4. Suppose now we start with the structure given in 
example 1, and we wish to change b(4;1) to 800 as was done 
in example 2. Suppose, however, that now we must maintain 
b(+;:1) and b. Under this restriction and following the re- 
sults of this section, we obtain the new billet structure 
B as 
385 240 180 70 55 
770455. 240 450 75 
Soom) 180120 60 
800 600 440 280 80 
Note that the first column of B still sums to 2340, the sum 
weenie ~irse column of B in exampie 2.1. Thus the drop in 
billets in b(4;1) has been distributed among the other 
focations. 
The new promotion matrix 1s 
poste». Ny 
7 oUGmn. 594 
Q = we AS, 


A446 eee al 
600 


Note again that only elements q(£j@and) p(l) change. ~Also 
mote that while in example 2.2 p(1) increased from .546 to 
molS, here p(t) has decreasiled to .517. The reason for this 
decrease is more subtle than was the reason fOr sei wLinGr ease 


mip t(l) im ¢xample 2.2. In the definition of p(1) as an 


40 





element of Q, it is seen that this number gives the fraction 
of personnel who enter a tour in grade 1 who subsequently 
EOD Lo@cmrnce tour in .crade Z, In the restructuring of the 
billet matrix B, we have increased b(3;1) from 300 to 385. 
Looking atmnemwe moce tLhatel(s; 1 )e-s ile whitle the other loca- 
tions have tour lengths of 3 years for personnel of grade l. 
Hence under the restructuring as given, more personnel com- 
plete a tour each year. Hence the "average'' number of pro- 
motions from grade 1 to grade 2 during a tour must decrease 
in order to maintain the same personnel flow into grade 2. 

In the mext section we senerafize theg§resultsi of 
this section to the case were all grades at a specified lo- 
cation ae scaled, reflecting a general increase or decrease 
Pomeeoe Gen Svzcmat the Specified location, 


3. scaling of Billets of all Grades at a Specified 
Location 


Suppose now that we wish to scale (up or down) the 
Population at yagspecified location, denoted &. In this casé 
we retain the proportional structure of the institution, but 
change the population at location 2 to reflect a general 
build-up or decay in operations? at the location. 

Henem vc amie Cm i (ee le othed) b Re eer rom Chis swe 
obtain b'(2;g) = b'(&3+)0(23g) = (1+a)b(2;+)0(23g) = (1+a)b(L3g) 
foreali= grades cae We also obtain b = b + ab(%3+). From the 


basic balance equations we obtain 
To a Cha eel eet oe Ko OL 
where 
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G 
Z 


K, = 124 D(&SK)R(ASK)w(k). 
We also obtain 
N, (k3a) = oe Sh 


vo NY, 
From (2.14) we obtain 


C, (k) + Ko" (k)o 
A(K30) = SEF BIER Gia” 


where 


Ko’ (k) = bB(Xsg)RE38) - 7, D(RsJIROVSI)IWG). 
We thus observe that the result is more complicated in 
nature than previous results. Note that q(k;a) is not lin- 
ear in a for any grade. Also note thatoq(k3;a)/da is not 
zero for any particular set of grades as it was in the 
previous work. This is expected, since now we have changed 
all grades at some location instead of a single grade as 


was done previously. We conclude this section with an 


example. 


Example 2.5. Suppose we begin with a system such as that in 
example 2.1. Suppose now that we wish to scale down location 
moyecetactor whien will bring b(4;1) to 800, the figure used 
memexanplieuw2z.26 oUppose alse that the,effect of the change 
1s to be only on the promotion scheme, leaving the rotation 
miccs wand withdrawal rates Constant. Thée billet structure 


becomes 
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500 240 180 70 Sa 
600 455 230 150 us 
500 240 180 120 60 
800 421 3509 196 56 


The resulting promotion matrix Q is 


. 354 ~ 546 
. 306 594 
Q = ~472 sees ; 
~491 2.09 
600° 

which we note is very close to the Q used in example 2.1. 
irs may oe explained by the fact that an across-the-board 
cut was made in personnel in the system. Due to the fact 
that all grades were affected in a similar fashion, there 
was not a significant change in the proportions of the 
grade structure. 
E. ABSORBING CHANGES IN ROTATION RATES THROUGH CHANGES IN 

THE MATRIX Q 

Suppose the rotation rate for a specific grade ata 
specific location is to be changed. Suppose also that 
enenec is) to béesno other change in rotation rates throughout 
the system. As was noted in sub-section D.1 above, in order 
PommMaintain the basic balance equations there will have to 
be changes in other factors. In this section we fix the 
Billet structure and the withdrawal rates and observe the 
Beecects on the matrix Q. 

Examination of the basic balance equations (2.4) and 


wero reveals that the-effects on the matrix Q of changes 
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in R with B fixed are symmetrical. Suppaee we are interested 
in the et fioeus or ta change in rotation rate for positions of 
moade feat tocation £5 1,.¢€. the quantity R(®3;¢g). Denoting 
meic=value of R(2iej}eatter the change by R'(L3¢e) = (ita)R(2 39), 
we note that after substitution into (2.4) and (2.5), the 
result shows the same mathematical relations ‘as were developed 
mm (2.11) and (2.12). 

Also in the same manner as was done above, simplifying 
assumption A2;4 can be made regarding the structure of the 
matrix Q, from which the changes in Q can be directly calcu- 
lated. Since there is nothing different about these calcu- 
lations once the vector Ny is known, we refer back to section 


ieee tor the method and compllete|) this section with an example. 


Example 2.6. Suppose that we wish to change the rotation 
rate R(4;1) to .234. We desire to leave the billet structure 
and withdrawal rates constant, absorbing the change in the 
promotion system. The resultant Q is that of example 2.2, 
aid since the product b(4;1)R(431) is the same in the two 
examples, this result confirms the similarity in the mathe- 


matical structure for the two situations. 


F, CHANGES WHICH PRESERVE THE MATRIX Q 
fm) Change in@assinghe Element of B 
Suppose the number of billets of grade g at location 
melse LOmbewechanecd. Jt 1s desired that the change not affect 
Pare Ppromotwon matrix QO. Since fixing Q implies that w and 


Ny are tixced, the Change must be balanced by a change in the 
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rotation structure R. Suppose only R(23;g) is to be allowed 
to change with the change in b(23g). 

With the above restrictions, we denote the change in 
D(Lsg)R(L3g) by b' (Ls g)R'(L3g) = (lta)b(L;g)R(L3g), obtain- 


ing functional expressions for iF and Ny as 


ro@= ro + wlg)b(4sge)RO38) , 


o(k) {r, + w(g)b(2sg)R(Qsg)a}-, 


N, (k3a) 
for k#g, and 


e(g) + b(X3g)R(L3g)a 
ry t D(Lsg)REsg)wlg)a 


Ns (g3@) 


From the restrictions, Ny (k; &) must be constant. Therefore 


p(k) - Nj (k;a)r 
a = - 
w(g)R(2;g)b (252) 


rN, CK) - rN, k3a) 
w(g)R(23g)b(258) 


fi 


= 0. 


Hence b(2;g)R(23;g) must be constant. Hence we are assured 
that the basic balance equations hold. 

Thus under the restriction of a fixed Q, as billets 
are increased, the rotation rate must decrease which means 
Maat the tour length 1s increased. Here*one can see the 
direct relation between the number of billets that can be 
supported at ajgiven location and the rotation rates or 
tour lengths. Conversely the above relations can easily be 
meveloped for changes in "supportable" billets with changes 


maerotautron rates. 
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Example 2.7. ees with the Orric ture of example 2.1., 
suppose we wish to change b(43;1) to 800 as was done in 
example 2.2, but this time we wish to compensate for the 
change by modifying the rotation rate R(43;1). Following 
the results of this section, it is readily determined that 


Res; T) 1s .475 for a new tour length of 2.105 years. 


2. Change in Rotation Rate R(2;g) to be Absorbed 
Throughout Grade g 


As we generalize from a single element change in 
R, discussed in section E, to a compensatory change through- 
out the system, the analogy of the rotation rate matrix, 
Or rotation structure, to the billet structure fails to 
hold. Although the mathematical structure carries through, 
the intuitive meaning of establishing rates which in some 
sense compensate for a specified change in the rate at some 
location for one grade is obscured. We specify the sense 
intended and develop relations which preserve the basic 
balance equations. 

Consider the effect on relation (2.5) if the rota- 
tion rate is independent of location, but remains a function 
Br orade,. Let R(*;k) represent this special type of rota- 


mom scheme, kK = 1,...,G. Then (2.5) becomes 


L 
R(*3k) 2, b(isk) = r)N, (Kk). 


But the summation on the left-hand-side is b(+;k). Writing 


Buc the right-hand-side of the above equation, we obtain 


Ilr 
R(*;k)b(+3k) = 42, bCisK)RCAsK) 
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il 


| L : 
or , R(*3k) = 52, bCisk)b(+;k) “RG3k), 


amd SR(*<kji ais seen to be a weighted average of the rotation 


rates. Further, using (2.8a) we obtain 


R(*3k) = 2, o(k;i)R(i;k). 
We are now able to interpret the constant Ky in equations 
(2.15) and (2.16). We see that for all locations & and 
grades g such that R(23;g) = R(*3g), the billet structure 
can be changed without affecting either r, (a) or the vector 
ere) - 

Suppose we start with a rotation scheme described 
Dy R(*3k), k=dl1,...,G. Suppose now that the rotation 
means 0) Lor a location £ and cyraide s91s to be varicd 
from R(*;g), independent of other considerations. We now 
develop the calculations which give a new weighted average 
rotation rate for positions of grade g at locations other 
than 2. This is done while at the same time ensuring that 
(2.5) holds. Denote the new rate R'(*;g). 

Since we desire that the resultant weighted average 


Mave the property that 


R'(*3g) 4,4. bCisg) + R(23g)b(k3g) 
be constant, we £ix R(*3;g) as the desired constant term. 


Bowe R*{*:2) can be wHitten as aefunction of R(2;9), giving 


R'(*5g5R(238)) = C, + KzR(X38), 
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where 


Ga 
i 


fee > 2) edits ye) (1 oO (E 1) ) 


and 


mn 
I 


3 = (gk) / C1 - O(g3k) ). 


It is essential at this point to check relation 
(2.4) to ensure that ty has not changed. Note that the 
above development is concerned with the changes in a spe- 
cific grade, denoted g. If ae Changes due to a change in 
R(%;¢) and R(*;¢g), then this change will affect all rela- 
tions (2.5) and not just the one in the variable g. Check- 
ing is a matter of straightforward substitution and it may 
readily be verified that ie does not change with the above 
scheme. Also since the calculations ensure that the product 
rN, (gs) remains constant and oe 1s constant, then N, (g) 1S 
constant. Silene + none of the other relations (2.5) are in- 
volved, then the vector Ny 1s constant with the prescribed 
changes in rotation rates. Hence the "promotion" matrix Q 
1s constant. 

Thus we have shown that rotation rates may be 
Changed so as to not affect other parameters of the model. 
One word of caution is due, however. As specific locations 
mor specmure 1OCation/grade elements arewsselected for the 
special" consideration regarding their rotation rates, the 
weighted average rotation rate assigned the remaining posi- 
tions may become unreasonable small (or large), implying 


unreasonable long (or short) tours. 
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An APL program was written which assumes the billet 
structure B and promotion matrix Q are given. The program 
calculates rotation rates based on the weighted average 
scheme described in this section. The program is inter- 
active in the sense that the user is asked to select any 
locations to be given "'special" consideration. Changes 
that result from the calculations are accomplished so as 
Eo not change the matrix Q or billet structure B. The pro- 


gram is given in an appendix. 


Example 2.8. Suppose now that we wish to change R(43;1) to 
.234, as we did in example 2.6. Suppose further that we 
wish to absorb the change by changing other rotation rates, 
holding other factors fixed. We obtain the new rotation 


mace matrix 


96 ee 55 OS, nS 


0. OE D Om are 
ies. ~ i i 1 1 
Boo ae ao oe 


which is observed to change only in the first column, i.e. 
the effect of the change is confined to the grade being 
changed. Calculation of Q using this R and the B and w of 
example 2 will confirm that Q remains the same under both 


motaction schemes. Thus the teffect ts felt only in R. 
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G. SUMMARY AND INTRODUCTION TO COMPARATIVE ANALYSIS OF 
ALTERNATIVES 


We now summarize the different approaches to rotation 
planning that have been discussed and generalize the results 
of the preceeding sections. The special cases used thus far 
will motivate the ensuing discussion. 

As we re-examine the special cases already discussed, 

a significant fact is evident. In the structure where two 
etry ibn = are jointly used in manpower planning, many re- 
strictions were required in order to obtain analytical re- 
sults. The underlying reason for most of these restrictions 
was the otherwise ensuing lack of uniqueness to the systems 
of equations generated. 

Whatever restrictions are placed on the model and what- 
ever scenario is thus selected for observation by the per- 
sonnel/manpower managers, there comes a time when one must 
ask which of several proposed alternatives is in some sense 
"best.''. In the usual case, there seems to be a hierarchy 
of factors considered in selecting the best of a given set 
of alternatives. For instance, the output of someone else's 
planning model has probably dictated the billet structure 
for the manpower planner. The most that the manpower planner 
can hope to accomplish with respect to altering the billet 
Structure is to point out significant problems associated 
with the imposed billets. For instance, there may be some 
hidden costs associated with one location. By pointing out 
these additional costs to his supervisors, the personnel 


manager may show that the operation in question is not as 
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profitable as the higher echelon planners thought. This 
approach may not always be necessary. Within his own purview 
there may be many alternative ways in which operations may 

be conducted. It is the manpower planner's job to select 

an acceptable alternative or set of alternatives. 

The sections that follow focus on the transfer cost in- 
curred due to a particular billet structure, rotation struc- 
ture and promotion scheme. It is to be remembered that the 
promotion scheme, identified with the matrix Q, may in fact 
be a retraining scheme or some other scheme involving the 
second attribute used in the rotation planning of the insti- 
tution. It is assumed that the measure of effectiveness 
selected to differentiate among the possible alternatives is 
the total transfer cost within the planning period. With 
these thoughts, we proceed into the evaluation of the alter- 


natives presented. 


H. A METHOD OF COMPARING TWO OR MORE PROPOSED ALTERNATIVES 
In this section we develop a method which enables us to 
compare alternative schemes that satisfy the basic balance 
equations (2.4) and (2.5). The method focuses on the mini- 
mization of the transfer costs between locations incurred 
during the period for which plans are being made, i.e. the 
average cost per enna in a steady-state model. We assume 
that the costs incurred in moving an individual of grade k 
meom location i to location 3, denoted c(1,j;k), are known. 
Now, let x(i,j;k) be the number of individuals of grade 


Peto be transferred from location i to location j during the 
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Pen vodweeletex(0i,);1) be the? number ae <ASane assigned to 
location j during the period. Here location 0 is meant to 
indicate an origination point outside the system, and all 
recruits are assumed to enter the system as members of grade 
1. The latter statement is a result of the assumption of a 
closed institution. We can then formulate the following 


problem, denoted P1, as a means of Comparing proposed 


alternatives: 
— * L : 
minimize z= fea SOs oc 2,5 1) 
L L G e . e . 
(P1) + i241 i Kz 4 S02 AES: CARRS) | (2.17a) 
subject to 
tf . | | 
52 SU STI) te eC R)) URIS } 
ae 
Ail eG de Oia (lk joe, el. y (27 eC) 
| 7= eee 
li 
i24 Dad sod eeh) Mees Kee |S Ligeees y 2 eka) 
Ke Jee. 4G 


and 
(eee) 


(COR) 8) Ey (ja = 2G 8D) aime on or 
Brere the xX(@e575k) are restricted to be non-negative. In the 
above, a(0;1) is seen to be the total number of recruits 
brought into the system. This is equivalent to the term 
previously defined as ive 
Note that the problem P1 is a transportation problem. 


Burther, mote that Pil is separable into G transportation 
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problems, since for any given grade k the constraints may 
be divided into mutually exclusive subsets of constraints, 
each containing only variables with third index k. We 

write out the kth subproblem, denoted Pi(k), for k greater 


eadh tle 


- L li 
Minimize 2z(k) = 24 fal Gian Westies 75k} (22 8ia) 


;2 
elit) ) 
| Subject to 


L 
j2y XC, j5k) = a(isk) i=1,...,L (2.180) 


and 


L i S © ® 
p21 *€4,53k) = r(jsk) j=1,...,L (2.18c) 
hiiere again the x(1,);K) are westricted to) be non-negative. 


For k=1 the above is modified because of the recruits, 


miving for P1(1): 
° L L « e ° ° 
min z(1) = jZq 42, C(i,js1)x(4,531) 


aei(1)) 
Sub ).e Gt le 
Ib 
aa ee = al yee 0 Le aul 
and 


L 
420 (1, eM Se 6 iL) a RS ie 


The problems P1(k) have a solution provided the right- 


hand sides, or constraint vectors, satisfy 


L ; Il 
p21 245k) = j21 r(j;k), k=2,...,G (2.19a) 
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and 


L 
= 421 1 Gage be (2.19b) 


© 
fo 
oo" 
Je 
Jed 
— 
| 


Equations (2.19) may be represented in vector notation as 


pe ee 
ayy EL ey 


21 pail (r.0,...,0), (2. Le) 


where 


L 
To = 52) x(0,j51) = a(051). 


Reference to section B quickly confirms that equation 
(2.19c) is the basic relation which led to equation (2.3) 
and subsequently to (2.5), one of the basic balance equations. 

Assuming now that the institution has the basic structure 
Piven in section B, we see that for every combination of B, 

Q and R satisfying the basic balance equations a feasible 
Solution to Pl exists; consequently, an optimal solution 
exists, provided only that costs are bounded. This optimal 
solution. will provide a minimum transfer cost for the stated 
B, Q and R. 

The optimal solution to Pl has EWiOlc LeMen Sammone £1St a 
Bost fiocure, a scalar representing the summations (2.17a). 
ma Ounce mec) CNet. BS) thesseL OLenuUMNDeLS On tmansters repre- 
sented by the xe (i ink). Note that the numbers x? (i,53k) 
do not serve to identify the specific individuals to be 
moved. The solution provides, in gross terms, the numbers 
of transfers of various types which would take place in a 
Sompletely controlled environment in-order to minimize the 


Sosts Of the rotation policy under the stated constraints. 
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We note that due to the separability, the vector x° (k) 
which is optimal in P1(k), when combined with all such vec- 
tors from the G subproblems, forms a vector x° which is 
optimal in Pl. 

One practical result of solving problem Pl can be shown 
by a dSeosqi nation of the number of basic variables in the 
solution to the problem. Pl has 2GL+1 constraints in L+L?G 
variables. Of the variables, at most 2GL will appear in the 
optimal solution at positive levels. In fact, due to the 
separability, at most G(2L-1)+1l basic variables will be in 
the solution at positive levels. Thus in an institution with 
10 locations and 9 grades, the optimal solution will have no 
more than 172 of the 910 possible variables that need be 
considered by the personnel managers. Thus a solution to 
the transportation problem greatly reduces the types of 
transfers which must be considered by the personnel manager. 

Suppose the manpower planner now wishes to investigate 
the effect of small changes in elements of the billet struc- 
ture or rotation rate matrix. We fone that under the re- 
strictions of section B the problemsP1(k) are related in a 
complex manner through the dependence of the elements a(i;k) 
on the requirements at location i for personnel of grades 
different from k and the matrix Q. That is, for any spe- 
cific availabilities a and requirements r that satisfy the 
basic balance equations, Pl will have a solution, and further 
Pil is completely separable into the subproblems P1(k). How- 
ever, for a change in either i eee or requirements 
Br ia na eeavetilen grade, other grades will be affected. 


De 





The dependence among elements of the right-hand-sides 
of problems Pi(k) leads to interesting problems. It also 
points out the desirability of many of the simplifying as- 
Sumptions made up to this point. Most assumptions made were 
necessary in order to obtain unique characterizations of the 
system. Without the assumptions, we are thus unable to ob- 
tain uniqueness. Without uniqueness, we are unable to spe- 
cify the effect of changes in requirements of a particular 
grade at a particular location on the other elements of the 
constraint vector to Pl. 

Given that we are able to uniquely specify the effects 
of such a change on the system, we wish to determine how 
these changes Shhce teat lomoo ti amecost ot Pie With the large 
number of parameters involved, we suggest some analytical 
wae Limb emcxpciacdmpmLor tOeeritccting the schanges, ihe 
development of metenods to effect this analytical effort is 


mie SUDJeECE of the next section. 


I. SENSITIVITY ANALYSIS OF THE TRANSPORTATION PROBLEM 
1. General Discussion 

In this section we develop methods which are useful 
meethne detection of “critical” elements among the right-hand- 
side of constraints (2.17b) through (2.17e). We begin by 
| considering the dual problem to problem Pl. We then develop 
meond -OrdGrecxpressions which yield the change?in the opti- 
mal cost of the transportation problem with changes in the 


controlled variable when feasibility. conditions are taken 
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into account. Specific expressions are dey elleped for the 
special cases cited in sections )D.1 and D.2 above. 

Consider first the dual problem to problem Pl, 
equations (2.17), denoted Dl. It is well known from duality 
theory that the value of the dual variables in the optimal 
solution to Dl provide a measure of the marginal change to 
z*, the optimal value of z, with respect to changes in the 
right hand sides of the constraints in the primal. Thus if 
u(i;k) is the dual variable associated with the constraint 
element af{2;k) in the primal problem, then the optimal 
value of u({i;k), denoted u? (i 53k), gives the change in z* 
with a unit change in a(i;k) provided the primal basis re- 


mains unchanged. Letting u(i;k) have this interpretation 


Lh 
v 


? 


mame MWettimsc v({5:k} Senic Gud. valiau lo asgociattea witn the 
& Jad) 


constraint element r(j;k), we can write the dual problem 


D1 as 


I 
maximize d = u({0;l)a(0;1) + 524 hey u(i3;k)a(i;k) 


Lee . 
ig ag VG aera) (2.20a) 


Subject to 


(Si eee (ys Ec (1.51) 1=0;5...,L (2. 20b) 
(a oe 
(lee ( Tek) c (1 ,) , Kein) , oak (2.20c) 
Heel. . 25h 
k=2, »6 


In this problem, the dual variables u and v are unrestricted. 
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Under the form of Q implied by assumption A234 we 
have 


ai eS) Peed Cie rl mem (w= tio) Ke2,..4.,G 


and 


aa ase a (te) q (1) 


Bon toe. be. oubstituting these relations into (2.20a), 
and recalling that assumption A2;3 implies that r(j;k) = 


b(j;k)R(j;k), we obtain, setting p(0) = 0, 


L G 
d = u(031)a(031) + 52) y2y vGsK)rG5k) 


(2.202) 
D 


G 
ea het 


MGs new soak yet © (i3;k-1)p(k-by 


ime sontimeal value for d is denoted d* and the values or the 
dual variables at d* are denoted u? (i;k) and v° (33k) fOr ga ie 
om jrand k, 

We note that Dl is separable, just as the primal 
problem Pl was separable. We write out the first and the 


kth subproblems 


Ik 
maximize d({1) = u(0;1)a(0;1) + 424 vai: Pag ly Rah) 
7 Ceci} 
D1(1) + 2, uG31)bG31I)RG31)q(1) 
SUD je CGtmmarommunad 1) + vi(¥3i) = c{(1,j3;1);52=0, sl 
ee, ore 


and 


a0 





L 
maximize d(k) = ity u(i;k){b(i;k)R(i;k)q(k) 


+ b(ij;k-1)R(ijk-1)p(k-1)} 
D1 (k) i222) 


+ 52) v(95k)b (5 5K) R (55K) 


Suvpeet. EO UCl skeet Vile; Kk =O C (tuo 5X) , eee 


BL foes 
| ES eee oe 

In the following subsection we use the above sub- 
problems in determining the effect on the optimal value of 
the kth subproblem of a change in billets of a selected 
grade. 

Peecilangves in a Single Element ot B 

Suppose the population of personnel of a specific 
prade, g, ait aq ampecatic lecation, 2, is increased or de- 
creased. Also suppose that there are no other changes in 
billets throughout the system. Thus the total population 
of the institution changes by the same amount. Suppose, 
as in section D.1 that rotation rates and withdrawal rates 
imemtixcd ss inus the buillet change 1s absorbed through a 
BnamGersli Ene matrix OQ. 

Recall from Theorem 2.3 that q(k) may be expressed 
as | 


1 
0 (k) 





(C, (k) = wigMR(2:elb(tsejal, k=1,..5,¢-1 


2 ty ay wOD0G) 
Bee aie “8) © ate) + Ritsa)bCsaya > *8 


C, (k) ile. Ge 


Nd 


a8, 





Thus we can observe the effect of a change in b(%3g) on 
q(k) by taking the derivative of the appropriate expres- 


Slons with respect to a. We thus obtain 


-j 
-w(g)R(2;g)b(L3;g)o(k) -, K Sse See 


oi : 
R(2sg)b(Xsg) ir, - «2, wj)eci)} 


dq(k3a) _ : 
30 Ce) * Rye) bs ga we 


0 »aogtl,...,G. 


With the exception of Di(g), and Di(gtl1), the subprob- 
Ronse given by {2.21 and (2:22) can *ee*expressed as functions 
of a by replacing q(k) with q(k;a) and p(k) with p(k;a). For 
Di(g) we obtain the optimal value of the objective function, 


d*(g), as 
L O 
d*(g) = ,2, u Cisg){bGsg)RGsg)q(gsa) + 


Diigee lL) Rie pte io) | 


L 
+ 52, vo (jsg)b(sg)RCsg) 


+ ab(2£3g)R(X3g){u°(R3g)q(gsa) + vo (23g)}, 


and for D(g+1), we obtain the optimal value of the obj ective 


munetion, d*{g+tl), ais 
I, on. 
Ce ee 2 taeee)yib{i eet yRGjetl)q(gtlsa) + 


b(i;g)R(i;g)p(g3a) } 
L O 
. f2i mes e tl) bi cme R(j5e+1) 


+ ab(£;g)R(Q3g)u (23 g+1)p(g3a). 


60 





We note immediately that for k less than g, d*(k) is 
linear in q(k;a) and p(k-1;a). These are in turn linear in 


a. Hence for all k less than g, d*(k) is linear ina. We 


obtain 
yd* (k) mae 
—zo = W(g)b(4sg) R38) a =Tky ; 
where 
L O 
Ay = w2y u (kj bp 5k Ria sk- 1) 
and 
: O 
By = 524 ode Dr ake) R (1s )e 


For the gth subproblem, we obtain. 
x 
on iB. Se ab(2;g)R(2sg)u°(2;g)} 2 Akio) 


+ b(<sg)R(L3g) uP (Lsg)a(gsa) + vo (23g) 


A_w(g) 


reese. 


Finally, for the g+lst subproblem, we obtain 


genet t) ye + ab(2;g)R(Rsg)u (X3g+1)} aa (esa) 


gel 
+ b(Msg)R(LQ3g)u (23 g+1) p(g;a). 


Hor all®subproblems indexed by k greater than g+1, there is 
no change in the objective function with a change in b(2%;g). 
Hence the above partial derivatives express the changes which 
mecurewitheaschangerin b(%;¢). 

Pole mlwovemeabculationssarescarried out for alll 


Peeotiencsmor the builet structure matrix B, the result will 
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be a display of the possible changes in transportation costs 
EO the various locations and grades. Abnormally high changes 
for some particular grade, location or grade/location com- 
binations would indicate that the particular community 


feserves cCloSeyscrutiny. 


J. RELAXATION OF SOME ASSUMPTIONS 
1. Relaxation of Assumption A2;3 (Stationarity) 

Suppose that we no longer require the requirements 
P(1;K;n) to be inde@pendent of n. Recall that in section 
B, assumption A2;3 was used to allow the basic balance 
equations to be expressed in the form of equations (2.4) 
and (2.5). We now develop a procedure for computationally 
allowing for the weak form of stationarity where the as- 
pip iwomeis that br(1;k:n) = r(i:kin-T(1i;k)). 

Recall that the above expression was required to 


mocain relation (2.3), 


Ib 
524 r(ij;n) = rj(@)N,, (255) 
and, under this assumption, (2.2) becomes 


r(ijkjn)w(k). an ee 


L eG 
To) = 42) aby 


Thus, for the system with "weak" stationarity but not station- 
arity, (2.2') and (2.3) become the system balance equations. 
TOmpmeccecdsenercetine QO(ijeto bewthe “promotion” 


Matrix during period n, from which N, (@) becomes the first 


row of the fundamental matrix during period n. Thus, under 
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assumption A2;4 on the structure of Q, a given set of require- 
ments r(i;k;n) and withdrawal rates w uniquely determines 
Olin). for those factors which satisfy (2.2') and (2.3), 
the transportation problem of section H may be defined and 
the transfer cost determined. 

The short-coming of the above is that the relation- 
ship between requirements r(i;k3;n) and billets b(i1;k) is ob- 


scured. This is best shown by an example. 


Example 2.9. Suppose we have an institution which has spe- 


cified the tour lengths T(i;k) by 


5 3 5 Z 2 


ON 
GN 
Se) 
DO 
DO 


and has identified the requirements for periods 1 through 6 


as 


200 160 120 80 40 
1 eed 6 ee ; 
300 240 180 120) 60 


M400 320 240 °&160 80! 


90 70 50 30 15 
210 150 110 70 35 

1 Ce!) | ; 
300 240 180 120 60 


*360 280 200 20 80 
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tl Ce 


ple 2) 


TR A), 


ee 


| 
wa“ | 
4 


300 
400 

90 
210 
300 
360 
110 
od 
300 


110 
190 
300 


380 


100 
200 


‘380 


90 
145 


20 


80 
160 
240 


280 


70 


240 


Zz 0 


90 
145 
240 
280 


70 
120 
180 
240 


60 
110 
180 
200 


40 

80 
120 
160 


50 
70 
120 
120 


40 
80 
120 
160 


30 
70 
120 
120 


20 
40 
60 
80 


80 


20 
40 
, and 
60 


80 


i> 
35 
60 
8 0 


The reader may wish to verify that r(i;k;n) = r(i;k;n-T(i;k)) 


for all locations 1, grades k and periods n. 


The reader may 


meso WIS Go Verity that r(::7/)°2 4r(;3;1), completing the 


eyclic. 


Suppose now that the withdrawal rates are given by 


wae (.1 6.3 2.2 


Bec t )i. 


Then we can calculate the Q(n) re- 


werred to Satisty the balance equations (2.2') and (2.3). 
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The result is given in Table 2.1. We can also calculate the 
recruits for each period, i.e. ro (n). We obtain the vector 
(660 604 658.5 616 648 614.5). 

However, the institution may wish to stabilize the 
eomoreten scheme, following the same promotion plan through- 
out the cycle. One way of determining such a "stabilized" 
promotion scheme would be to determine the promotion matrix 
Q based on the billets within the institution and the given 
rotation scheme. This will be demonstrated through the use 


of another example. 


Example 2.10. For an institution with the structure used in 
IG eS ae 


— 


example 2.9, the relation r(i;k) = b(1;k) gives 


the billet structure of the institution as 


300 240 180 70 35 
600 455 230 150 US 
300 240 180 120 60 


‘1140 600 440 280 80 


as may be Beachy, verified. Thus the basic structure of this 
example is identical with that in example 2.1 and we have a 
"stabilized" promotion matrix from that example. Using this 
approach, we expect some disparities between availabilities 
and requirements in the short run, but would hope that they 
mmerace souteover the lone run. Table 2.2 gives the excess 
mt availabilities Over requirements for the six periods and 


for all location/pay grade combinations within each period. 
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.340 .560 
.300 .400 
Q(1) i= .466 .333 
: .500 .200 
600 
.371 = .529 
.314 .386 
Cw) = .476 .324 
AW Omen A 
600 
528) eee S72 
£295 .415 
0 (3) = .472 .328 
.500 .200 
600 
.384 .526 
.321 .389 
Q(4) = .476 .324 
A / Gus 2 4 
609 
- 7.325 .575 
7 7 AOS 
Q(5) = Gil GES 
.500 .200 
600 
.373 .527 
.316 8.384 
Q(6) = ue a wale 
EG 1, Baas 
600 


irple 2.1; ihe Promotion Rates for Periods 1 to 6 for 
Example 2.9. 
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b(C 2 5 IL) ee 


nA 52) 


E32) % 


Requirements 


Over 


nemexaGesscsSmon Wva1babili tives 


able 2.2. 


tou rPervodse 1 through 6 for Example 2.10. 


67 





Pieatablo=Zs2 we iawe deiined E to be the array of thevexcesses, 


where element e(i;k3;n) is a(i;k;n) - r(i;k;n). We observe 
that 

li G 

i2y Kz1 e(i;k;n) = “ro (1) 


for all periods n. Referring to section B, we see from the 
discussion leading up to equation (2.3) that this is to be 
expected. 

We nNote@that the elements of Table 292 Ware rela- 
tively small compared with the requirements for each period. 
To measure the differences between the scheme used in example 
2.9, where the promotion matrix was Galeutaced. £orscach 
period and the scheme used in this example, we let Eg be the 
Si idy OLMeXcesisicls an example 2.9 and LO be the array of 


Seo soc Smo Ven minim Lable 27. lhenewe Galleuliaate the  per- 


centage differences between the two schemes, denoted P, from 


C19 (i3k3n) - €g (13;k5n) 

p(i;k3;n) = a OO aa aa x 100. 
PpleeZ  seprvesmune array Pa FromyP wesmote that the largest 
meLreentaye crror 1s approximately 3.8, occurring at location 
mecouefTadces during period 6. Wihus aie least for thiseex- 
ample, ie net effect of all the changes in promotion 
Schemes amounts to an increase in ability to meet require- 
Menmts OL Only 3.8%. 


eee laxXaimon of Assumption AZ;1§(Single Q) 





SUDpOsecailow tilat we assiime that there exists a matrix 


Q. providing the transformation on past requirements at 
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1.36 1.09 0.34 i / S Z. 16 

P(::1) = I ASKS 1.09 0.34 1.78 2, Lom 
Ges 1.36 1.09 Oy 34 i. 78 ZO 
UES: 1.09 0.34 1.78 246 
Ley S ie 47 Jey ke! ie 1.2 254 
P(332) 273 ie Ow d Kees ae ie. 1s Ue es 
pre 1 ier 12 4a 1.09 3 Zsa 
7/3 7 i 4 1.12 1.91 
eS S 2.0% a eS. 0.94 216 
Dia = (eS ss leg 9k Lage By 0.96 BG 
Ar 75 5 2.08 56 0.96 2.16 
2 255 1.92 F.56 0.96 2.16 
5204 ge) 1.49 ZeLG 2.54 
Pern ie Su 4 BESS) 1.68 Zaew7 ies | 
om 3.04 250) 1.49 WS Poa 
fe) 4 28 5 1.59 Zale 1.91 
LAS 2.50 0.89 Zou Ze. 1.6 
P(::5) = 2.86 205 0.69 Ones 2. ZAG 
NE 2.86 Teg 0.80 Za 2 2.16 
2.86 iG 4 0.80 Zo 2 2.16 
1.94 1.38 O18 5) tele) peo 
npeee) = 1.94 56 0221 2.98 2.54 
Med 1.94 iy 4 4 O23 2.90 2.54 
1.94 io 5 0.29 3.08 reo! 

mole 2.5. ihe Percentage Difference Between the Excesscs 


Ce wAvaglabilities#Over Requirements Under the 
POCO s Twin | Cm. 10 sieskxample 2.9. 
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location i giving current availabilities. This would seem 
to be a "natural" approach to the problem. Differing tour 
lengths at different locations would seem to result in dif- 
ferent promotion opportunities at different locations. We 


now have 

a(i) = r(ijQ. SL 1 ease 
or for the basic balance equations 

le ik th 

52, 2G) = 42, r€i)Q; = 52, rj) - (rQ,05--.,0), 
and 


L G 


i234 k= 1 r(i;k)w(i;k) = r 


Be 
meereswiigk) 1s agfunction of location since Q. 1s. 

Thus the set of promotion matrices Q.> i= eee 
and requirements provides the availabilities and recruits. 
Unfortunately, much of the analytical work derived in the 
preceeding 1s lost to us. For instance,sgiven the require- 
ments r and the withdrawal rates w (now a matrix) we are 
Mnablie to obtain @ unique set of Q. from the above two 
relations. 

JemGe velo pmunGuditbecreme appradenes=to the problem 
with location-dependent matrices Q.. First, suppose for 
each location 1 we are given a matrix Be such thait B00 ~ Q. ; 
fooecllSigand a Common promotion matrix Q. Recall that the 
Pymbol © implies)clement by element multiplication. Under 
assumption AZ2;4 we have, letting sie be tire (},k)th 


a oe 
Poment Of the matrix E , 
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L 
«2, r(isk){1 - e*(1,1)q(1,1)} 


a 
1! 


and 


0 = a Ma sk)i{i - Be alae) 


- r(i;k-1)e"(k-1,k)q(k-1,k), 


the latter statement holding for k=2,...,G. These equations 
can be iteratively solved for the elements q(m,k) of the 
matrix O- 

The generalization seems to have some merit, but the 
determination of the matrices se may be difficult. Using 
this in the model for the previous work may lead to results 
that are heavily dependent on particular eo used. Since the 
elements of the ee do not seem to be directly measureable, 
we do not pursue this course. 

Another approach to the problem with location depen- 
dent promotion matrices Q. 1s somewhat more general. The 


Straightforward generalization of (2.2) gives 


L G 
ro (m) = 524 Key rani ule (cia sk) ) Ww (a7 ke (2.24) 


where the modification expresses the withdrawal rate as a 
function of location as well as grade. 
We also have, under the assumption of weak station- 


petty, that 
(2.250) 


gO) a SoCs t jn) = 1 r(i;m;n)q, (m;1)}, 


Shovel ‘ 


(Al 





L G | 
O= .2, trGiskjn) - 2, r(ism3n)q, (m;k)}, (2.26) 


for k=2,...,G. Note that in the latter two relations, the 
elements of the promotion matrix Q.; are indicated by the 
subscript. | 

These equations may be manipulated in much the same 
manner as was used in sections D, E and F above. We provide 
a sample of the calculations involved. Suppose that the 
requirements at some location and grade are to be changed. 
Suppose also that no other requirements are to be atfected 
and that only the promotion matrix for the specified location 
is to be changed. This would be the case, for instance, 
where each location in the institution is autonomous insofar 
Pe rromotions are concerned. Of course, another inteérpreéeta- 
tion would be that if the rotation rate for a particular 
location/grade were decreased (i.e., the tour length becomes 
larger), then the promotion rate per tour should increase. 

SpecitiGadilly let us) assume that the requirements tat 
mocatliOnmum.tor persionnel of grade g are to be changed. Let 
the new requirements be denoted r'(%;g) and be given by 
fate OO) = (ita)r(23s). finger the assumption A2;4 of no 
demotions and at most one promotion per tour, we can calcu- 
Mette Ciesettect of thisWehange on Qo» SayMee BULIT eyclikere Q. 


remain fixed. 


Let dy (k3k) =—les B(k))q, (k;k) for all the diagonal 


elements of Qy- We obtain 
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B (k) aes, c- L 


ar (% ;g)w(%5g) k 
r(L3K)qy (k3k) ” 
ad, (g5g+1) 
Pts) = (Teayq, (ese) ’ 


and 


B (k) 


I 
© 
Me 


| Clee Be, SG 


Bhis approach providesean interesting alternative to 
Perce approach slelected for use in the bulk of the’ chapter. 
'Computationally, the magnitude of the required storage is 
seen to be much larger here. Rather than one promotion 
Matrix of size GxG, we now have L matrices of the same size, 
one for each location. Also, instead of a 1xG vector of 
withdrawal rates, w, the withdrawal rates are now specified 
by a LxG matrix. | 

Byereneile@e Of One approach over the other seems? to 


be application dependent. 


K. EXTENSIONS AND CONCLUSIONS 
1. Summary 

In the foregoing, we have approached a significant 
Manpower problem - that of the interactive structure between 
billets, rotation rates and promotions. Under assumptions 
Piiciewere required Jn ordér to obtain Uniqueness and thus 
dial ytitcal tractability, we have obtained a structure which 
us versatile and useful in providing insight into the man- 


power system of a large, multi-lJocation institution. 





The versatility of the structure has been demon- 
strated through several sections which gave some of the dif- 
ferent approaches possible to effect a change within the 
system. A convenient means of evaluating the alternatives 
presented by the different approaches was given in the form 
of a transportation problem format. The manner in which 
various elements affect the cost of the transportation prob- 
lem was discussed. A method of detecting "critical" elements 
of the transportation problem was developed. 

Pee UXTeHsS bONS 

There are several possible extensions to the above 
development. We observe a few of the more significant. 

a. Multiple Time Periods 

it may well be that current decisions in the 
rotation scheme will affect future rotations. Each person 
rotated during this period must rotate at the end of the 
assigned tour, Either out of the system or to a different 
location. Thus we should be concerned not only with current 
costs, but should also consider future costs associated with 
errentsecxpenditures. This)may be done through a multiple 
meme periodstransportation problems It may be desirable to 
Brscount the future costs incurred asja result of current 
MecisSions. 

SUCH SanMapproach might Mivert the cffects of 
sub-optimization, i.e. making short-run "optimal" decisions 


which are costly in the long-run. 
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b. Restricted Flows 

It seems reasonable that many of the possible 
transfers would not be desirable. For instance, in the U.S. 
Navy one reason for rotations is the maintenance of equitable 
periods of sea and shore duty among the personnel within a 
Skill group. In this case a transfer from Stein Diego-Sea to 
Norfolk-Sea duty might not be allowed. 

In the preceeding discussion this problem was 
implicitly resolved by assuming that high enough costs c(i,j;k) 
could be assigned to these transfers to drive them out of 
the optimal solution. This of course assumes that there is 
some set of the x(i,j;k) which is feasible and doesn't include 
any flows restricted by the special considerations. Other- 
mise, the optimal cost is inflated due to the necessity of 
remaining feasible and the associated high cost of elements 
which are restricted. 

Hopefully, the use of sensitivity analysis would 
Lav cmLOMmITdiCatcwany Vidlationm of SUcCh @ restricted arc. 

If in an actual implementation the Peeeniett ons on flows do 
cause problems, it may be worthwhile to consider them ex- 
plicitly. This may mean that one will have to forego the 
enice’™ structure of the transportation problem. 

Geb ans 1 Ciimmcad Sc 

Although much of the above work was involved with 
Eee cmangime Of thc personnel istructure of the Institution 
gn various ways, the results were obtained from a steady- 


State model. Thus the model shows the magnitude of changes 


fis 








in going from one structure to noche The short-run ef- 
fect of a change in the system was essentially ignored. 

It would be of interest to develop relations, 
for instance, to show how to "optimally" change the promo- 
tion scheme in a sequential nature to get from steady-state 
under one billet structure to steady-state under another 
billet structure. Included would be the sequential changes 
in the billet structure which would be a consequence of the 
sequentially changing promotion scheme. 

Since most large systems probably never acquire 
steady-state operations, such developments would be helpful 
in determing how to get from the current personnel structure 
of the institution to a desired structure. 

d. Stochastic Considerations 

In the preceeding, it was assumed that the 
Hiller structure and rotation Paecse wWetemolvcnemoOreda lterna— 
tively that availabilities and requirements were given. It 
was also assumed that the system followed the established 
meructuremexactly in@all) casies!. Umnethe usual case, only 
estimates of system parameters are available and the system 
will fluctuate in a non-deterministic nature about the Spee 
eetted rotation scheme, 

For anstanececachmetencnte= aqi;k;n) of the 
muatlabi lity, vector might be considered to be the result 
mies@arch Ofsthe personnel-files (giving the number of per- 
ponmnel ofserade k at location iv’who are to become available 


during period n) amended by such considerations as the 
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fraction which will be promoted to k+tl prior to or during 
period n, the fraction who will leave the system and thus 
not be available for reassignment and the fraction who will 
Demtransterredmpetween thewtlmewoet) thessearen and the period 
n. Estimates of these fractions might be obtained through 
analysis of historical data, examination of legal contracts, 
promotion policies (current and projected) and/or managerial 
intutition. 

Similarly, requirements are projected based on 
the results of a search of files, amended by similar con- 
Siderations. Many factors which are less predictable than 
those given above influence the actual availabilities and 
naa renen te of the system. One of the more significant is in 
mle ailvcumepwea tour. A tour iS Saia to béseinterrupted if the 
individual concerned rotates in a period other than that pe- 
riod in which the maximum tour length for his present posi- 
tion is completed. There are many reasons for interrupted 
tours. For instance, an individual may require hospitaliza- 
tion and/or convalescence for a lengthy period of time. This 
may force the institution to assign the position to another 
individual. For humanitarian reasons such as death or 
serious illness in a family, an individual may desire relo- 
cation prior to the completion of his present tour. An in- 
dividual eligible for retirement may decide to retire prior 
Mem previolsly expresscd iantentions’. An employee may be ac- 
cused of commission of a crime and thus be unavailable while 


awaiting trial. Someone may die. 
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The above describes some situations in which 
tours may be interrupted due to events in the personal lives 
of employees. Institutional decisions may also cause tours 
to be interrupted. It may be determined that operations at 
some location are to be curtailed and personnel reassigned. 
Some positions may be temporary in nature; with the attain- 
ment of the objectives of a particular temporary position, 
the individual in that position may be reassigned. Simi- 
larly, operations may be initiated or expanded. Temporary 
positions may be created. All of these decisions lead to 
availabilities and/or requirements which might not be in- 
cluded in the rotation scheme. The early availabilities | 
are directly a cause of interrupted tours. The unforeseen 
fequirements may cause interrupted tours throughout the 
mictieution aS individuals are selected to fill the 
requirements. | 

Actual availabilities and requirements ance an 
mectedshby factors other than interrupted tours. Personnel 
may transfer between skill groups, possibly maintaining 
their pay grade. See, for instance, Hayne [21]. These 
Pransicrs sdremtcrmcam lateral transfers. Wlso, the institu- 
BinitemaveUSe@tne concept of®al'igapped” billet, 1.e€. a posi- 
mioiewhtcheexists On paper, but need’not be filled”at all 
times. Obviously, gapped billets provide some flexibility 
mimene Tequivements vector. Also, the transiént problem 
discussed above might leave some dislocations within the 


eysocen when moving from one steady-state structure to another. 
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Thus, invalid estimates of see ROO Unpredvet= 
able fluctuations and changes in institutional policy can 
lead to dislocations in the rotation scheme. Such disloca- 
tions may cause severe variations in the cost of the rota- 
tion scheme from that given as a solution to problem Pl. 

Particular schemes for relating availabilities, 
requirements and the relations between the two while at the 
Same time accounting for fluctuations might be the subject 
of significant statistical study. Upon development of the 
Statistical relations, revised estimates (and/or distribu- 
tional descriptions) could be used in the formulation of a 
stochastic transportation problem. Since during the period 
of consideration all availabilities and requirements become 
known, the structure seems to fit that of a stochastic pro- 
gram with recourse, see for instance Walkup and Wets [36] 
and Williams [38]. 

The stochastic nature of the elements of the 
problem would be applicable if the model were to be used as 
a forecasting tool. The present use as a planning model can 
be accomplished quite well without the stochastic considerations. 

e. Personnel Planning 

| iie=prececdinmop chapter hase been concerned with 
manpower planning. That is, the idealized system was struc- 
tured to obtain results in terms of the over-all flows, etc. 
On the other hand, personnel planning 1s usually construed 
to imply concern with the ‘unex AR Which to manipulate 
Proolunencoources in Otder to best achieve the "goals" of 
fie institution. 
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By slightly restructuring the above model, it 
could be useful in personnel planning. For eatin Comat 
the availabilities and requirements are taken from a search 
of the personnel files, then the transportation problem would 
give the gross numbers of transfers between locations to 
give the least total transportation costs for the period. 
Since the optimal solution must be feasible, the solution 
to the transportation problem also ensures that if personnel 
are distributed under the scheme in a deterministic system, 
then all requirements will be met. 

Such a structure could also include stochastic 
elements; thus the comments of the preceeding section apply 
here also. 

3. Conclusions 

We. see that the structure developed herein has rami- 
fications throughout the manpower and personnel planning of 
adelarge multi-location institution. We have been especially 
concerned with the interaction between billets, EOCAt1om 
rates and promotion schemes as they aftect the rotation 
problem. We have formalized the structure among these ele- 
ments and provided a means to compare the effects of various 
elanges on the resultant transportation cost of the institution. 

Our cffort 1s seen as a beginning, the comments a-e 
miesub-2scction 2 above indicating some of the directions that 
io tibcmuncwmr taken ss (he pursuit of better and better undcr- 
standing of the complex interactions between components of 


aiespersonnel! system of large institutions si of increasing 
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interest. Diseconomies caused by failure to understand these 
interactions are of concern to both institutional directors 


and their constituents, both public and private. 
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eel bae oe ROONNE tmokieCl LON.PROCESS WITH @RANSFER COST 
CONSIDERATIONS 
A. INTRODUCTION 

In this chapter we develop a model of the personnel 
selection process within a large institution, where selec- 
tions are made to determine which individual transfers are 
to be made. In the previous chapter we determined the gross 
numbers of transfers of personnel of grade k from location 
mecombocation j aS the x(i,j;k) in the transportation yprob- 
lem. Here we are concerned with which one of the a(i;k) 
individuals of grade k who are available toc be transferred 
from location i will be sent to location j to satisfy a 
meoguumemomn amonrmtie r{(q;Kk}opetully, ¢he sclections 
Will be made so that the total number of transfers during 
the planning period will approximate the x (i,j 6k) from the 
peevious model. 

The process to be modelled has as its principal concern 
iow 1NclinGero nde dUd ladeicdamep baCeCnecnumtouttliagSingle given 
demand. Here qualified might mean satisfactory performance 
Mim OGereduisivewpOSIULONS, Of 1t might*mean attainment of a 
Epeclilespayecrade, Oouettamioht mean qust finding a ‘body" 
momid!lgaeslot. ) Ihe degree of qualification will not be 
occ med tur the r- 

The. general approach taken iis’ to observe thé behavior of 
Mn idealized system. The idealized system to be studied is 


PS tbictedtsin Scope and application, but lends itself to 
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providing insight into the structure of the problem. Later 
sections generalize the results somewhat. 

The first system we study we term the "Ehrenfest Decision 
Model."' Here we draw on the results of the well-known 
Ehrenfest Urn model of heat flow between isolated bodies 
(for instance, see Feller [9]}. The known Pesidts are 
adapted into a decision-making framework. In the special 
case of the Ehrenfest Urn problem, one becomes concerned 
with the choice set over which decisions are to be made. 
Results are obtained in the form of probability distributions 
over the choice set. From this, distributional results de- 
scribing the transfer among locations are obtained. 

We develop some approximations which are helpful a that 
they provide insight into significant parameters of the 
system. 

The Ehrenfest Urn model per se is usually considered as 
consisting of two isolated bodies. The Ehrenfest Decision 
Model used here is developed for the case of two locations, 
mictecencraivzed to three locations. )— Further, the choice 
Set 1S. generalized to include cases where sequential deci- 
Sions are allowed. 

A more general type of generalization is a description 
of the time dependence of the Mis caiticas. This leads 


directly to the use of diffusion approximations. 
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Pee ineebaotG MODEL += THE TWO LOCATION PROBLEM 
ie Hin endest Decision Model 

Consider a closed firm or institution. The concept 
of a closed institution and the various reasons for trans- 
fers between the locations in a multi-location institution 
were discussed in section A of chapter II above. Suppose 
now that the institution has only two locations where opera- 
tions are conducted. 

As in chapter II, we assume that personnel are de- 
scribed by two attributes. Here we use current location 
and qualification as the two attributes of concern. Thus it 
may be that we are restricted to personnel of a single skill 
category/pay grade combination. We may require also that 
Preto temileGualltiacatlons be Cyvidencedgwithin the rescriccrea 
group of personnel to be considered in the selection process. 

We assume the existence of a central decision making 
organization, commonly called an asstgnment center or per- 
sonnel distribution center, whose function is to determine 
which individual to assign in response to a requirement or 
demand. 

We call the two locations R (for red) and B (for 
Pack )wemNemcall the personnel distributiongcenter M* (for 
manager). Demands made to M* for personnel originate at R 
Siebe they also include a specification offthe qualifica- 
miOnmereduiredeto pertormethe duties of the position. In 
mrceconte:.t Of chapter II, demands occur due to forthcoming 
Bemretions Of Cours, tours being interrupted, or simply 


creation of new positions. 
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Costs are assumed to be incurred whenever a replace- 
Mcitui@ce ma pemsionsselected for transter nm response fo 
the demand, is currently in the location other than the one 
originating the demand. In this sense, a physical separa- 
“ser, Aes the locations is assumed to exist. The institution 
is assumed to be liable for the expenses incurred in moving 
the replacement, his family and their possessions from his 
current location to the site of his next position. 

Unecwagdcnand 15° received at the personnel? distribu- 
tion center M*, it may take considerable time and effort to 
locate a suitable replacement for the individual currently 
in the position. The assignment center must find an indi- 
vidual who is now, or soon will be, available for reassign- 
ment and who meets the specified qualifications. 

The assessment of qualification may include an ap- 
praisal of each individual's past assignments and perform- 
DiCCumEENCEDUnpOSc: ©. the rotation plan must be considered. 
i, for inistance, one of the farm's goalsefor the rotiation 
program is to provide exposure to varying aspects of the 
firm's operations, the assignment center would not want to 
assign someone to a position he has previously held. Other 
factors might be important. The promotional opportunity 
provadedeby a particular posithon’ may correlate with the 
promotional potential of individuals assigned that position. 
aiateus yea high exposure’ position may only be piven to 
SNewot descleect group within the broad classijficetion being 


Pansiacred. 








Some theoretical consideration of the time to find 
a suitable replacement when searching a list of potentially 
(altoid personne! 2S givem ain Chapter IV; 

Upon determination of a replacement, the person cur- 
rently in the position for which the demand was generated is 
assumed to become available for reassignment. Thus his per- 
sonnel file is added to the group which is searched by the 
UiStriputton center upon reccipt ofga demand. Thus)every 
exit from the group available for reassignment is matched 
by an entry into the group. 

Here we ignore any interactions between groups 
aVidtlavleeGrOreheasSHenment. t.cC. eroups oL different yskill 


types or pay grades. Thus the size of the population avail- 


HWemecomstructesa simplegmodel wretlectine the structure 
developed above. Suppose the records of those personnel 
available for reassignment are considered to be balls ina 
"decision" urn M*, where M denotes the fixed total number 
Sooballsene the Urn sinew balis am M* are painted’ red’ or 
mhacn, dependinces on the cumnment physical location of the 
andividuail whose fiile is represented by the ball. 

Let Rin) denote the number of red balls in urn M* 
PO lOWIMG@ the nth transfer to occur in the sysitem. Simi- 
lanly B(n) denotes the number of black balls in M* at the 
SS poOtiumiin tinCnmelcarniy, R(n) + Bin) = M for all n. 

Demands may be considered to be the outcome of in- 


q@ependent Bernovlli trials. This is reasonable in the sense 
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that if there are Np positions of a certain type at location 
R, and Np positions of the same type at location B, then the 
fraction of demands originating from R should approximate 

Np/Np + Ny). We also consider demands to occur independently, 
and let r be the probability that a random demand originates 


from location R, where 


We also let b = 1-r, i.e. the probability that the demand 
originates from location B. 


Graphically, the model may be represented as Figure 


eee Ce a iG 


ae >) SC al 
transfers 





Pores wee lhe Iwo Location Ehrentest Decision Model. 


In the next section, we obtain results on the steady- 


Peace mols ityedis tribution of the composition of urn M*, 
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2. Distributional Results Under a Simple Decision Rule 


Suppose the personnel distribution center M* simply 
finds the first qualified person and effects a Bayeten. In 
the simplest case where all personnel whose files are repre- 
sented by balls in M* are equally qualified for all demands, 
this would imply that the policy is to take a random draw 
PeOm UlrieMioand €ffect, aatranster. 

In more complex situations, there may be a lengthy 
search to find qualified personnel. In the latter case, the 
search would terminate upon the location of the first person 
qualified for the position from which the demand originated. 

In thise section, thessteady-states probability dis- 
PripiiitonwoLeki(ijels) cosbe determined. Ihe value of R{n) ait 
any time n is called the state of the system. Also to be 
determined 1s the steady-state probability that an arbitrary 
move 1smad COSt Move, lew represents a transfer “between 
the two locations. 

Let R({*) denote the number of red balls in urn M* 
Prteheanlwabvltrarllyslargesnumber of transfers have occurred, 
Thus R(?) is a random variable describing the number of red 
balls in urn M* once the system is in steady-state. Thus we 
Wish to ‘find the probability distribution of R(~). Note that 
Hetermination of the probability distyybution of R(®) also 
PVCs ticmdeterninatvon of thesprobability distribution 


oo B(o): 
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Theorem 3.1. Under the above conditions, R(*) is distributed 


Binomially with parameters M and r, i.e. R(~e)~Bin(M,r) 


Before proceeding with the proof of theorem 3.1, it 
is useful to develop some notation. Let Pr {A = a } be the 
probability that the discrete random variable A assumes the 
value a. Let Pr { A = a | B= b } be the conditional prob- 
ability that the discrete random variable A assumes the 
value a given that the discrete random variable B assumes 
the value b. 


Proof: Consider the conditional probability statements: 


Pr { R(n) = j[R@-1)=j-1} = SU ,, 

Pr { R(m) = j[R(n-1)=j } =dr+ Dhar) (3.0 
and 

Pr { R(n} = j|R(m-1}=jrl } = i= (1-r), 


mnere the (1rst Cquation state's that if R{n-1) =) 3-1, 

R(n) = j 1f and only if a demand arose at R (an occurrence 
with probability r) and a black ball was Chosen from urn M* 
Pe Ccenlrnrence with probabiimty [-(j-1)7M). Similar anter- 
pretations apply to Chew OIereevGse<presslonses Note that 
changes CL MOnLemelan once ula tmare excluced, 


Noweeletne(z) be three probability peéneratinge function, 


R(n) wy 3. 
Ceo) = ket 2 lige 524 limi = i fs. C552) 
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_ 


where E{A} is the expected value of the random variable A. 
At steady-state, Pr{R(n)=j} = Pr{R(n-1)=j}. Let p(j) de- 
note Pr{R(o)=j}. Assuming steady-state, we —_. the 
condition on Rin-t) in equations’ (3.1), obtaining 


a(z3e) = 3Zy 27 (a rp(j-1) + Z#2ca-r) pen) 


+ { if. + it (1-r)} p(s) } ’ 


or 


doz 52) er 


ae och Zz 8 (z3%) 


i 
© 


The solution to the above differential equation subject to 
8(1;°)=1 is 
8(z;”) = (1-rtrz)". 

Hence R(~) is distributed binomially with parameters M and 
c: 

ine nex tegucSt1Ol WGmmcOncCe ml 1oucnatroL Ene prob= 
ability of an arbitrary move being one between locations. 
Let c(n|j) be the PLoOo bie midi wemeloOcat lon. branS- 
BoreOccurs on the nth itransfer piven thateR(n-)-j). «We 


obtain 


e(nfj) = (1-1) + Fer 


r+ d(1-2r). 


At steady-state the unconditioning in straightforward and 
People Mettingse(*) be the probability that an arbitrary 
fieve US Ome inter-location move when the system is in stéady 


state, 
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(oy a cele yr) 


Thus if C is the cost of transfer between location R and B, 
the expected cost of satisfying an arbitrary demand when the 
system 1S in steady-state is 2r(1-r)C. 

In the next section, we model a more .complex system, 
where two or more qualified personnel are selected prior to 
a decision being made as to which one to transfer. Distri- 
butional results are obtained for the steady-state case. 

Swe variable Length Selection bist, (5) 

Suppose now that the personnel distribution center 
is required to find a specified number of qualified person- 
Ife OnLome tO ettecting a transter, Denote by & the number 
specified. Ignoring momentarily the problems associated 
with the search time to find the L specified people, we 
analyze the Ehrenfest Decision Model in a manner similar to 
the above. in this section, however, the state-dependent 
Markovian properties of the system are useful in the analysis. 

he ttranster policy@specified for this section is 
EPiateiieany on the L selected personnel are currently located 
at the location from which the demand originated, then one 
of those is transferred to the new position and the remain- 
ie eeepetsounel! tiles! dre retaurned togM*am Hencegan inter— 
location transfer occurs only when all L selectees are 
Gurremtly in the location other than the one from which the 
demand originated. 

bee tomonovrols that Some restrictions cxisit on L. 


Pomeiactoncomiemie were. Jargecr than M, the size of the 
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population in M*, then the list would never be completed. 


The relationship between L and M is now investigated. 


Theorem 3.2. (a) If the selection list has length less 
than or equal to ie then the system in steady-state is 
absorbed in an absorbing class with M-2L+3 states. 

(b) If the selection list has length greater 
than or equal to ie then the system steady-state is 


absorbed in one of 2L-M-1 absorbing states. 


Hoos, (a) Le< 1+> implies tit #- ie Male oe Bie Rin) < bat 
ieneRintl) > R({n). That?is, each@selection of?L balls?when 
R(n) < L-1 must have at least one black ball. Hence under 
mic tEransier policy stated above, mor R{n) in the specified 
mon Comore bell wrill@eyer be transSerred te B. Bence 
R(n+1) can never be less than R(n). 

Pik ects ieethen Bin} belies by varguiemt simi 
Pi mecomeme mDOVer (nt lig-s b(n). Hence@R(net1) = R(n)- 

Furthermore, in each case there is some positive 


probability that the inequality will be strict. For 


Rin) < bL-1, 
M- Jj 
am i) 
Pre (RG) > Gy GMS a 51 oS IB ab “hee 
| Gi 


A . ; Ag. A! 
where (p) denotes the combinatorial, 1.€.(p) = 50 Gear 


Por eveercatucr than or icqualyBe For R(n) > M-L41, we obtain 
J 
ae: (7,) 
Pr { R(m+1) < R(n}|R()=j,j > M-L+1} = — (1-r). 


M 
(1,) 


WZ 





Hence R(~) may take on values only between L-1 and M-Ltl. 
Puncthormorege it bsereadily seen that tor L < R(n)*< M-L, 
there is positive probability that R(nt+1) = R(n), that 
R(n+1) = R(n)t+i1 and that R(nt+1) = R(n)-1. Hence ail values 
from L-1 to Morea Agee Ene Gancepot R(@) 2 Hemce thes sys- 
tem is absorbed in an absorbing class with M-Lt+1-(L-2) or 
Paro SbaAtes ; 

(b) L> 1+5 implies L-1 > M-L+1. I£ R(n) > M-L+1 
for any n, then the argument above gives us that R(n+1) < 
Fle Aliso, if R(m) < L-1 we have that R(nt1) >= R{m). 
Hence for M-L+1 < R(n) < L-1 we have that R(nt1) = R{m). 
Hence each state from M-L+1 to L-1 is an absorbing state 


and we have L-1-(M-L) or 2L-M-1 absorbing states. 


Pxaipecmoa ieee LCteM=11, L=558 The one?step transition matrix 


ms represented by 


R(nt+1) 
R(n) Grete ABS home Oe oe Lonel 
0 xXx 
1 b, aoe 
Z x xX 
5 Dow x 
4 Xx. x 
5) (> xe x 
6 pene O. 
i ne 2 
8 xX 
9 me OX 
10 x xX 
Lil xa xX 


Pie TOmuLAwea ore presciu positives pmovabilities. To trace an 


evolution towards steady-state, suppose R(n) = 2. Then R(ntl) 
oi take om vValucs Zor 3. I# it asisum@s the value 2, the 
nex#® transfer 1s under the same probabilistic structure. If 


we 





Riise chem Rin) > 2 tor all m > n+l, by the above 
theorem. The system continues operation until 4 < R{m) < 7 
for some m. (Obviously in this case, the absorbing class is 
entered whenever R(n) = 4 for the first time.) Once in the 
absorbing class, a probability distribution can be deter- 
mined which will describe the state of the system at an 
arbitrary time in the future (i.e. after many transitions 
have occurred). We do this. 

For the rest of this chapter, L is considered to be 
less than i + > This is the case of interest in that it 
lends itself to the determination of a steady-state prob- 
ability distribution. The other case discussed above is 
obviously dependent on initial conditions and hence is of 
lamiGed interesc. 

We write out the conditional probability statements 


comparable to equations (3.1) for the case of variable length 


selection lists: 


Pr { R(n)=j|R(n-1)=j-1 } = 


J eaig { R(n)=j3 |R(n-1)=7} = ] AL (1-r) = maha Tr, (sey) 


and 


j+1 
a) 


Ei RGN 3) ales thas } = —_ fil-7), 


Cre, 


ie MOO SSE Tt wnat @thcr transitions 1s excluded. 
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From the proof of theorem 3.2, we know that states 
Ujt,...,6-2 and sitates’ MFL+2,...,M alre transient. To deter- 


mine the steady-state distribution of M*, we solve 
T= 1 P, 


where 7 is a 1x(M+1) vector with elements 7(j) = Pr{R(~)=j} 
BOG js0.1 oe M, and P is#the one step transition matrix 


with elements ao Pr{R(n)=j|R(n-1)=i}. We have that 


m(O0) = e¢* = wW(L-2) = W(M-L+2) = *** = m(M) = O, and 


M-L+1 
j2L-1 ™(j) hes 


io Ome mrOL Mom TL) S| reCUrsiVely, first Obtaining an ex- 


pression for each t(j) in terms of m(L-1) as 
” 6 ) 
m(Ltk) ie) Jy-—_————— n(L-1), k=0,...,M-2L+1. 


Mi Toe : 
fettine Ado = JH, ath Or andeecy =, a 7 eh ero 


relation is written as 


k+] 
| ms A (k) ; 
1 (L+k) = (5) O16) TCL ie we (Sieccs) 
M-L+] Mea Ltt 
But we also have that ,2) 4 m(j) = m(L-1) + 52 1(L+j) 


fieiecmmsubostitutine (53:4) anto the above, 


M-2L+1 
- Pan P Atk) 
meee ( + = (-) ate \o 





Hence, 


Wea) Sor ay , j=0,...,M-2L41. 


M-Z2Lt+1 
k+1 k 


1 + Lu (ts) 


A 
C ( 


The preceeding relation is not in a computationally efficient 


form. Letting D(j) = A(j)/C(j) we see 


M-L+1] 


D(O) = —+—, (3. 6a) 


(7) 





and 


ee 
it | 
ee) 


D(1) 


HI 





DD) 


Paegeneral we can write 


M-L+1-3 
: Li : : 
DOF = i Lae, Do ee) 0 eee Ce 
(bt) 
OT 
D(j) = (3.6b) 





mommaeseoe M(M-Z2L-9) = D4), 1=0 , ee, | - L, 


where [ ] indicates "greatest integer less than or equal" the 





waLuc indicated. 


Proce by Straichitforward isubstiytution into (3.6), we have 


p(o) = —Guszbttj! 0! 


(Magi 1)' Lt? 
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and 


L!(M-2L+1)! 


eS eso (M=be1)? 


D(M-2L). 


By evaluating A(M-2L+1) and C(M-2L+1), we see that D(M-2L+1) 


= 1]. Hence 


eae = Ue] je 


D(M-2L) = tM-2ne1)2 LP 


= D(0). 
Hence the lemma is true for j3=0. Assume the lemma is true 


Mal, 





Zor Some j within the allowable range, i.e. j=0,..., I. 


We show that the lemma holds for j+1. 


(M-L-j)!(j+1)! | 
W205" (ej+tyy PO) 


D(j+1) 


D(M-2L-j). 





But from (3.6b) 


eee (Lt jal.) (MaZL-7)! ~ ee 


By transposing, we see that 


D(M-2L-(j+1)) = D(j+1). 


Hence: theslemma is ishown by induction. 

The %signiticance of the lemma is in the reduction of the 
number of coefficients D(j) which must be calculated in order 
PomedGtermine the steady state probabilities from equation 
5) 

Calicut tons to détermine. the 47(L+j) ts! are streightfor- 


Perteonee tie D0(,)'s are computed. A FORTRAN program was 


oF 





to do the computations. A sample of the results for M=20, 
1 06 |p ae Oe 2 tie and. tor re".1,.2,...,.5 are fiven in 
Pabbes nl and S22. 

We note that for L=1 the values in Tables 3.1 and 3.2 
are binomial which agrees with the theoretical results of 
Scetmon Baz above.» in fable 3.1, the effect on the Steady- 
state distribution of urn M* resulting from an increase in 
L is shown. Such a table would be of interest for a stable 
system with a known r. If the system under study were to 
Contemplate changes in the personnel Structure, then Table 
Seczewould Show for a given L the effect of a change injr on 
the steady-state distribution of M*. 

The next section develops some approximations to the 
calcuiations of this section. These approximations are 


used to determine the effect of changing L on the steady- 


state probability of incurring a cost move. 


4. Approximations to the Steady-State Probabilities of 
the Variable List Length Process 


In this section we develop some approximations to the 
steady-state probabilities developed in the last section. 
Comparison values are computed for the same conditions 
undem the procedure of the last sectiom andy under the 
approximation. 

feeLeusemuch fesis than M, as seems reasonable tin Pmany 
Case's, then one might treat the model as sclection with re- 
Gecceicimmmundereselectrom with replacement, the one step 


PELensbemone probabilities are baised on powers of the fractional 


vie 
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R(*) = 


R (2) 


Mable 3.2. 
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M = 20 L = 8 continued 
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The steady-state probabilities) for the two location 
variable Length yselection iist@problem, M = 20. 
Hommine specitfaled Wength of scileictiton Tiist, L, the 
Probabilities are given as functions of r, the 
probability that an arbitrary demand arises at 

leca tion he 


104 


a 





composition of urn M* rather than on the combinatorials 
encountered in the last section. Equations (3.3) are ap- 


proximated by 


| L 
Pr{R(n)=j[R(n-1)=j-1} = (G2), (3. 3a) 
— e L 
Pr{R(n)=j|R(n-1)=j } =1- (ii) (l-r) - (“7 ) Le 
en! 
and Pr{R(n)=j|R(n-1)=j+1} = (Ar) (1-r) 


Letting p(j3;n) = Pr{R(n)=j}, we obtain 


ee) 16) | Ghee R ele) a a) (1-r)p(j+1;n-1) 


P(j3n) 


+ (1 - (i) (1-r) - (Mi r)p(j;n-1). 


The above equation is representative of a system of 
equations which are indexed on the states j and time periods 
n. These equations are quite messy. <A simple deterministic 
version, suggested’ by Gaver, may be obtained” by letting p(t) 
denote the number of red balls in M* at time t. Time is 
considered to be measured in units of transitions. Then 
we }can fexpresis "the time ratte of change of p(t) in terms of 
oe Clmee rite traction of red *baliisein urneM= at time®t will 
increase if L black balls are chosen in response to a demand 
PEO andewill decreasco if L red balls gare chosen in re- 
sponse to a demand from B. Otherwisc, the fraction of red 
Daliceewmd Sthust the fraction of black balls, will remain the 


Same. Jhus the time nate of change of p(t) may be cxpressed 


by 








Speed = (Sree wy": (ect) (1-r). Gap, 


\- 


At steady-state, denotalp = p(), oe = 0, which im- 


plies, letting £ = p/M denote the expected value of the 


fraction of red balls in urn M* at steady-state, 


Ges: & eis). 


Thus, thesexpected value of the fraction of red balls in urn 


M* is 


le. 


f = : (3.8) 
WE. (ery 27h 


r 


Thus the quantity Mf should approximate the Saal number 
feeuedmbal isp im M* at gsiteady-state. We mote that for L=1, 
the approximation is exact. 

Tapilcmer smod Voss Valllcsato: the=expeceeuw number Of red 
balls in M* at steady-state as determined by the approxima- 
Elon MEyepmccedoedsey thesgexact value Calculated from the 
probabilities determined in eC rion 3 above. The approxima- 
tion is seen to be closer to the true value for small L and 
hOmevalucseoLt rr near 25. Thésresult yaor smald L is a result 
Seechesassumption that L isemuch less than M. In Table 3.3, 
mhe Urnesaae Meis 20, tand thle values of®L near l0eane cer- 
tainly not small compared than M. Also, as r approaches .5, 
miemCmpeetCa 1 racthon of red bal Ms’ 1m M* approaches”; as? can 


pesscen fromplables 3.1 and 3.2. This coincides with the 


moctUnhtia( 5.5) Lor r= .5S. 
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ME ME ME ME ME 

i, = a 2-4 4 6 6 g g 10 10 
2 sie a 6.76 6.67 7.97 7.92 9.02 8.99 10 10 
eecaicsf9 7 89° 7.75 8.70) 860 9637 9.352 10 10 

4 7.66 7.32 8.50 8.28 9.08 8.94 9.54 9.49 10 10 


5 BR io Or mo O72 eed Se oe 7D LO S10 

6 8.64 8.19 9.13 8.85 9.46 9.29 9.75 9.66 10 10 

fi S.C Omebe ooo, o2 90 9.57 9559 9.80 “9271 10 10 

8 O17. 8.65.59 47 9245-5 9.608 VO ey 9-785 9.75 10 10 

9 COMET, CMO OU MO m2 5 9, Oo. oS eee.ieo 9. 77 LOG 

10 9.56 8.90 9.75 9.31 9.86 9.58 9.94 9.80 10 10 
Table 3.3. The Expected Number of Red Balls in Urn M* at Steady- 


State, Followed by the Approximation as Given in 
(3.8). The Figures are for a System with M = 20. 


Te sd lich Grse Between Locations sim tne Variable List Length 


Dele emlol t nocess 

In section 3 above, distributional results describing 
tiems Usady-Statee composition Of Urn M* wereWalichieved. Given 
those results, the movements within the system are of in- 
Cres te 

We first determine the steady-state probabilities 
Eiatyan arbitrary move Us agtransfer between locations, 1.e. 
is a cost move. We then develop two approximations to the 


desired probability statements. 


Oe 





Letting c(n|j;L) = Pr {nth move is a cost move given 


that R(n-1)=j and the selection list is set to be L}, we 


obtain 
M-j j 
Cy) 
c(aljst) = Ber + Maen, 
OP GE 


Now letting c(L) denote the steady-state probability that an 
arbitrary move 1s a cost move when the selection list has 


length L, we can remove the condition on R(n-1) obtaining 


ee Mk k 
; Goliad! (1) 

e(L) = HY aes. uf (1-r) Te Kj} eo toe os 
=e ee 1) 


Once again, we can resort to a simple deterministic 
approximation. Let c'(t;L) be the approximate cross-transfer 
Rae .l,c. sche rate otfetlow between the two Jocations. Then, 
using the arguments leading up to equation (3.7), we have 

a By Mane) L o(t) L= : 
e'(t;L) = Cae. | G: ; (i-. 


At steady-state, p = o0(*) and oe = 0. We have 


2e(1-r 


(3.10) 
Me Gly eee 


©) (CUNT oame vat 


[r : 


otc titi Or Le=5 i, the above@equation gives the résult 
MEA Ned In Section 3. 
Die waperoximation €‘(1,) can be further approximated. 


Wesurite c'(1,) in the following sjecial way: 
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ies a sg) 


eb) = ae tC. 
ae _ Lert “-1) + leet ) 


We now observe the following relations. Let x = (rib aay, 
Then pith = ce + le Then r =@(1 + ox ae Thus as L grows 


1. 
large, r approaches exp(2x) or x approaches &n(r*). Sub- 
stituting this limiting result into the above equation for 


c'(L), we obtain 


eo) & r(i-r) : 
2 i ZL 
ob 1 (1 -, SMT) : 
L 1 
But the second term in the denominator can also be approxi- 
mated by an exponential, resulting in the second approxima- 


fom, cgenoted c''(L); 


aveiaye nee @eiul 


The last approximation provides some useful insight 
iinomene Dehavior On stheasystemeas 1 HS vVarieds= Ihe cross- 
transfer rate is seen to decrease by a factor of 2 GE Le 
L increases. Thus, for L=2 approximately half as many inter- 
Neca bon MOWVeSmare 1equared asi were required for Lai. Main- 
taining the movement rate for L=l as a base, we obtain that 
Pomc om tietemire approximate lye.eas many intcer-location 


moves, and for L=4, there are approximately 1/8 as many 


moves. 
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a 


Thus L is seen to have a profound influence on the 
cross-transfer rates of the system. Chapter IV develops the 
difficulties encountered at the personnel distribution cen- 
ter when L is increased. We observe here that the cost 
Sane obtained by decreasing the cross-transfer rate does 
introduces difficulties in other parts’ of the system. The 
magnitude of these difficulties must be weighted against 
the savings to be expected by increasing L. 

Table 3.4 compares the two approximations c'(L) and 
c"(L) with the exact cross-transfer rate c(L). Examination 
of the table shows that the approximations are not good for _ 
L greater than 6. However, the purpose of the approximations 
was to gain intuition about the influence of the length of 
Bice sa lc cmon LS ion themerass-transten rate. Ome sees 
immediately that the approximations under-estimate the in- 
fluence of L on the transfers. Thus the result above that 
the number of transfersi, or the cross-transter ratte, de- 


(L-1) understates the influence of L. 


creases by 2. 
Piet hemiLo<tSeCet non Wwe cenchall ze theseirentese 
Decision Model discussed in this section to the case of 


three locations. 


See lLieeciReNe con DECISION MODEL WITH THREE AND*MORE BOCATIONS 
imeeeccneral Duiscussi on 
The generalization of the two location Khrenfest 
Decision Model to three and more locations is of interest. 


Thé large institutions which are of concern to us throughout 


EEO 





c(L) c'(L) c"(L) 


L=ti 420 420 458 
Z Pea fag he, n229 
3 099 = i sito 
4 042 056 A017) 
9 .017 .028 O29 
6 006 .014 014 
i .002 007 .007 


Table 3.4. Comparison of Two Approximations to the Cross- 
Transfer Rate with the Exact Value. The entries 
are for a system with M=20 and r= .3. 


the thesis would normally have personnel stationed at many 
locations. The effect of various peroonnel pos Teles =o tine 
expected costs to be incurred in the execution of those 
policies is the subject of investigation. 

In this section we are principally concerned with 
Piece cate Gt one@hne length of thet sicleetion list som Ene scosts 
incurred. We begin by looking at the trans fers @vnich soccur 
when the system responds to a demand by assigning the first 
mUalidied person that it 1s able to locate. Thiiseis the 
i teease.s We then imtroduce I as a parameter am the 
description of the selection process. We find that in the 
generalization preferences of various types of transfers | 
MUScmmemeS tiated. for auspeciigc preference, we develop the 
mathematical relations for the one-step transition matrix 


and then develop steady-state approximations. 





We then investigate the costs incurred for the 
variable L case. We are unable to obtain analytical results 
in this case, but numerical results can be obtained. We 
conclude the section with a discussion of the complexities 
introduced by the generalization. The scope of the added 
complexities is formulated, giving specific results on the 
Size of the one-step transition matrix and the numbers of 
calculations which are required to obtain the elements of 
the matrix. 

2. The Three Location Problem with L=1 

Suppose now that the institution under study has 
three locations at which personnel are situated. Except 
for this modification, the system is assumed to behave in 
the same manner as that described in section B.2 above. 

Here we have arbitrary demands originating from one 
of three locations, R (red), B (black) and G (green). The 
demands are processed by the personnel selection center, 

Or manager, again denoted M*. Upon receipt of a demand, M* 
Searches the records of the available personnel in the spe- 
cified grade until a qualified replacement is found, where- 
upon that individual is transferred to satisfy the demand. 
Upon completion of this process, the individual in the posi- 
tion from which the demand originated is made available to 
M* for reassignment. 

To define the composition of the records at M* fol- 
Powine Ghee nth transfer, we specify the vector (R(n),B(n),G(n)). 


(it welchentsmonetmce vector give tiféenumber of records at M* 


pi. 





which represent personnel physically located at R, B and G, 
respectively. Again letting M be the total number of 
records at M*, i.e. the total number of records of person- 
nel who are available for transfer, we have M = R(n) + B(n) 
ot ia). Recall that, in the two location problem, specifica- 
tion of R(n) was sufficient to describe the composition of 
M*. In the three location problem we must specify the cur- 
rent values of any two of the elements of the vector 
(R(n),B(n),G(n)). Choosing R(n) and B(n), we define M(n) = 
(R(n),B(n)) as the "composition vector." 

LC Gaemioees atid g be the probabilities that an arbi- 
trary demand originates from location R, B and G respectively. 
Thus we have that r + b + g=1. A graphical representation 


Peo tic VrOmIciH Or ChiS¥SCCLION IS)g1Ven in figure 3.2% 





aos cemands 


See Dhysicak 
trans leis 


} i sentati > Pocat von 
Figure 3.2. A Graphical Representation of the Three Loca 
Ehrenfest Decision Model. 


BS 





We have the one-step transition probabilities 


pr{M(n)= (j,k) |M(n-1)=(3,k)} = dr + Ky + Mork, 





M M 
Pr{M(n)=(j,k) |M(m-1)=(j-1,k+1)} = ae 
Pr{M(n)=(j,k)|M(n-1)=(j-1,k)} = U-U-k, 
PriM(n)4(j sk) |M@el = (7*1, (k-1)} = It (Gua 


Il 

= 

| 
<i 
ro 

| 
Me 

Is 


Pr{M(n)=(j ,k) |M(m-1)=(j,k-1)} 





PriM(nje Gk) Mas etki te ie ; 
and Pr{M(n)=(j ,k) |M(n-1)=(j ,k+1)} = a 


Graphically, the possible states for urn M*, where 
M ="5, may be represented as shown in figure 3.3. In this 
aoc wet hnomemiackened duke weet Sent the possible feasible 
states of the system. Mathematically, the feasible states 


i emt he mlon-navative integer solutions) ta the equation 
yjeteke + 2 =8M. 


Note that the j-coordinate represents the number of red balls 
iimurh M eenom k= COOrdildmwc represents Ene number of black 
balls in urn M* and the &-coordinate represents the number 
ei@creen Dialls in t@hewurn. 

Note that the planar surface given by j + k +2 =M 
La mcSstricted tO thcmpas tive Orthant by Vimtue of the 
nonmenegative requir@ément, 1.e. that there cannot be a 
negative number of balls in the urn. 

EStmBii ks) —eRriM(mj=a(j,k)}. Then from the system 


Sie equiactomomio. 12), Wwe have 
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Figure 3.3. Graphical Representation of the Possibie States 
of Urn M*, when M = 5. The blackened dots 
represent the possible states. 


p(j.ksn) = Cdr + Sp + ME *y yp Gj ksn-1) (3.13) 


+ Med yp (j-1,k+15n-1) + dbp (j+1, k-15n-1) 


Mo) kt la +a cs | 
+ en) Erp (j-1,k3n-1) “ +7 -ep(j+1,k;n-1) 





a Meda p (j sk-13n-1) + Shep (j,k+1jn-1). 


M 
Define 
| | (3.14) 
M M : : 
poe ee = R(n) B(n) B: lads ene? < 





Substituting (3.13) into (3.14) and assuming steady-state 


conditions exist, we have 


@(Z45253%) = 0(24.2)) 


, TZ4 + bz, ued - ies 
(1- (rz, +bz5+g) )M , a 9Z4 


C2) as 


90(24525) 
re) 


The solution to this partial differential equation subject 


EOeg (l,l) 1 ys 


0(2Z,52Z5) - (124 v bz. + gyi (3. 1S) 
From (3.15) it is immediately recognized that the 
composition of urn M* obeys a multinomial distribution with 
parameters M, r'and b where g=1-r-b. 
As in section B.2 above, the next question of concern 
is that of the probability of an arbitrary move being one 
between locations. Here we let c(n|j,k) = Pr {inter-location 


move occurs on the nth transfer given that M(n-1)=(j,k)}. 


We then obtain 


: B apr + kb ft (M- jak 
c(n|j,k) =1 - 227 Me 


from which, by unconditioning on M(n-1) and assuming steady- 


state, we obtain 


(aa 0 cee = be 
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We note that for the two center problem, WOEtlINo Belay an 
the above fashion, we obtain c(~) = 1 - r* - b*, from which 
we readily obtain c(@) = 2r(l-r), which was the result ob- 
tained in section B.2. 

Comparing (3.12) with (3.1) we note the increase in 
complexity that results with the addition of ‘another loca- 
tion. The degree of increase in complexity of the mathe- 
matical description of the system with increases in numbers 
of locations is discussed in section 6 below. 

As inithe cassie of the two center problem, the next 
extension is to the case of varying the selection list 
lee tic nilseis Gie SUbJECct Of the nextesection: 


3. The Three Location Problem with Variable Length 
Sel e.G i aOnmeeers © 


Im this section, we tstudy the three locaittion EBhrentest 
Decision Model under the conditions of a variable length se- 
Mcectron liste, fhe length of the list is again denoted L: 

We see that the analysis becomes much more complicated than 
the analysis of the two location problem with variable L. 

Pi etiam te=Ocatt One (imeem orcater than two loeca- 
tions) problem, some preferential ordering must be included 
Mimo amc: Cc ihe Lc ephObl Clim lo CXdMimesthe need ror 
Pit SwOLdecrimes Suppose 42 demands ormpanatecs fatelocation R, 
Further suppose that in the selection of candidates to fill 
the demand, a selection Jist of ‘length L ¥#s' required. Upon 
ays@arch of the) available personnel, the list of candidates 


may include personnel currently located at B and G, but none 





located at R. A decision must then be made to send one of 
the candidates at B or one of the candidates at G. 

If the decision is based on the minimization of 
the, transfer cost incurred in satisfying this one demand, 
then the required preference ordering is easily determined. 
Suppose that in the institution we are studying the transfer 


costs between centers have the following relationship: 


Cop 7? Cag 7? Cpe 


where Cop is the cost of transferring one individual from 
HOcietonek to locations B. =) themother costs@ifave similargin- 
terpretations. Under the assumption that we are to minimize 
the cost of each individual transfer, the ordering of pre- 
ferential transfer may be represented by table 3.5. In 


table 3.5, lower numbers indicate preferred types of 


mransters. 


LOocaeion, Ou (Candiaqatre 


R 1 4 g 
= 
ZB A 1 2 
=D 
 G 3 2 1 


Wable 3.5. The Preferentilal Ordering for Transfers in the 
Three hoca mon variap ie? List lene th Problem. 
Lower numbers indicate preferred transfers. 
For instance, a demand originating at location R would 
preterably bewmet by a candidate currently located at R. 
Powmlemomednatante trom R be included in the L personnel 


Be receved them a Gandidate from G (if found) would be sent 
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to R to satisfy the demand. Only in ight case’that all 
selected personnel are currently at B would a person from 
B be sent to satisfy the demand. Similar interpretations 
apply to the remaining elements of the matrix. 

| Graphically, the system may be represented by 
figure 3.4. The numbers on the dotted arcs indicate the 


pacrecmence ordaecringwon the arcs. 





—— > demands 


——-—--— physical 
transfers 


Figure 3.4. A Graphical Representation of the Three Location 
Ehrenfest Decision Model with Variable Length 
Selection List. Numbers on’ arcs represent the 
preference ordering on transfers. 


We can write out the conditional probability state- 


Men esmeompanapie to (5.5). Wevhave 








Pr {M(n)=(j ,k) |M(n-1)=(j,k)} =] - == - Ls - ty 
ee) <aaw 
‘wie 
Pr{M(n)= (j,k) |M(n-1)=(j-1,k+1)} = “44, 
ae 
(nae) On 
Pr{M(n)=(j,k)|M(n-1)=(j-1,kK)} = —ye— [ 1- gee |r, 
=a) 1 
Ce i 
er 
Pr{M(n)=(j,k) |M(n-1)=(j+1,k-1)} = +, (3.16) 
Gs) 
rs eee 
Pr{M(n)=(j,k)|M(n-1)=(j,k-1) } = —H—(1- CH Jp, 
Cy) L )/ 
jr 
PriM(n}=(.k)|M(n-1)=(+1,k) } = SR , 
(gn) 
and L 
mie) (4) 
Pr{M(n)=(j,k)|M(n-1)=(j,k+1) } = —~e[{1- =e |g 
jt+k+] 
aw) des 


To interpret the above equations, let us examine the 
FILSGe second and third equations above. “lhe first equatoions 
simply states that the probability of remaining in the same 
State 1s one minus the probability that no reds were selected 
times the probability that the demand was from red minus 
SUC tmeiise WL thmiesiMiiktn intcrpmeration.  IThhnis cquation 
describes ali type 1 moves, i.e. the most preferred. 

Diese conameduecnon tn (S.16)etypifties a type.3, or 


least preferred, move. Here the move described is from B to 
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Rew | feBequatron states that this?move occurs only iff all 
candidates selected are currently at B and the demand 
originated at R. 

The third equation in (3.16) is slightly more com- 
plicated in nature. It 1s an example of a type 2 move. 
The move described is from G to R. This move occurs only 
when a demand originates at R, there are no reds in the list 
of candidates and that of the balls selected, not all are 


black, Note that in the third equation, the statement 


+a 


M 
H) 





wells Suse ties probability that nonesof the j-1 réd balls in 
M* were selected. Also the statement 
G) 


an.) 


1 - 


mellis us§ the probability that not all®the balil’s’ selected (of 
which we know no reds to be selected) are black. Thus the 
denominator in the last term reveals that the population 
bemigmesamnpilede 1s jonly of8sazeeMa(j-1), singe the population 
available to sample is noun to Bewnon-rcd. 

We are unable to obtain explicit analytical results 
POumeie WiljsmnrOMeT wa Te, since anspection Shows that from 
Anearbiltrary starting point aSemany as iseven states can be 
reacned, one for each equation which applié¢s. Howéver, in 


aepanrerctiiar application the one-step transition matrix P 


a 





can be generated from (3.16). Then the t(j)'s can be ob- 
tained numerically. 
In the next section we develop approximations in a 


manner similar to that used in the two location problem 


above. 
4. Approximations to the Steady-State Probabilities 
sf thegMaria blew Liste lero ther recess ewathimibre ec 
Locations 
We again develop some approximations to the steady- 
state probabilities given by (3.16). As was done in section 


B.4, we suppose that L is much smaller than M. This allows 
us to approximate the above development by using the con- 
Sot wun crecltion wath replacement equattons (5.16) are 


then approximated by 


Pr{M(n)=(j,k) |M(m-1)=(j,k)} = 1 - iy te a Ky D (Ky g, 











L 
Pr{M(n)= (j,k) [M(n-1)=(j-1,k+1)} = GA) or, 
L 
fe YL 1) = ey pe) | MiGs = (7-1, k)7 2 (=) G Cr Satie Je 
; ; jel z 
Pr{M(n)= (j,k) [M(m-1)=(j+1,k-D} = GAA) b , 
| | eye ile — 
eMC i) IM (Co = Oe ke (Gea) of B 
ey 
Pr{M(n)=(j,k){M(n-1)=(j+1,k)} = (4) g , 
a , 
aya Pr{M(n)= (j,k) [M(n-1)=(j,k+1)} = (at ———— ) “(eget g 


BACs 





The above equations are seen to be in powers of L 
rather than in combinatorials aaatoladiae L. At this time, we 
could uncondition on the M(n-1). However, we saw in section 
B.4 that the result was not very helpful. 

We proceed to a deterministic version. Let p(t) 
dencte the number of red balls in urn M* at time t. Also 
let B(t) denote the number of black balls in urn M* at time 


t. Let time be measured in units of transitions, Then 





do(t) . (8(2)) 


+ (ROY (2 - (SEB) ) = 
: (oft). (b+g) , 


and 








Au 
Lo») 
laa 
c 
Nm 
ll 
P aa 
Le) 
Fe 
cr 
Nee! 
” 
ar 
+ 
ons, 


| at L 
sa) (- (w8tB) ) 
Ear (o(t)*s(t)! ( (ergot ) Je 


Letting f(t) denote the fraction of balls in urn M* 





which are red, i.e. f(t) = o(t)/M, and h(t) denote the frac- 
tion which are black, h(t) = B(t)/M, we can rewrite the above 


equations in the fractional form. We obtain 


th 


meet) 


sell = (12 - £(t) )¢ r - £(t)*(1-7), 


and 
y tht) . (1 h(t) dp - h(t) rr + F(t)’ 


at 
fj) Ween we 


TS 





The latter equations may be obtained using the following 
argument. We observe from Table 3.3 that location R is 

the "least preferred" locativn from which or to which trans- 
fers are to be effected. That is, the preference ordering 
is established such that the higher values are placed on the 
interaction of R with the other locations. 

Thus, transfers into location R occur only if no reds 
are selected in response to a demand at R. This event occurs 
with a probability of Glee Geh ioe. Also transfers from R oc- 
cur only if all the selected balls are red and the demand 
arose at a location other than red, an event which occurs 
with probability £(t)'(1-r). 

Similarly transfers into location B occur only if 
no blacks are selected in response to a demand at B, an 
event which occurs with probability Coeanwee Sb. Transfers 
from B are of two types. From Table 3.3 we see that a trans- 
Per tromesd tamke1s sere least preferred offall types of trans- 
fers. This occurs in response to a demand at R, and only 
if all the selected balls are black, i.e. with probability 
h(t)?r. Titemo thei. PC mG tT anS te hat Gomme bed Smiter CSpOnse 
to a demand at G. If the selections made are all either 
red or gacall (i.e., no green balls selected) and there is 
at least one black ball selected, then a transfer from B to 
G occurs. This event occurs with probability (£(t)+h(t))H¢ 
- Clete 

At =steady-state, the derivatives are zero and the 


iiemoOnsmoicanondonger time-dependent. Letting £ = lim f(t) 
t->oo 
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and h = lim h(t), we can then write 
t+ 


ames. Boas Ss 
and : (eS 


L 


(1-h)"b - h’r - (£thyg + thg = 0, 


The first equation has the solution given by (3.8). This 

is expected since the preference ordering of Table 3.3 
established location R as the location with the least inter- 
action with the other locations. Thus the other locations 
could be considered as being one entity and the three loca- 
tion problem reduces to a two location problem insofar as 
transfers to and from R are concerned. 

Unfortunately, we are unable to obtain a closed form 
Pim Stee tO S.3) for the sceondtdquation im (3717 ji. 
In the next section, we investigate the transfer costs in- 
curred in the three location problem. 


2 iranstters Between Locaitionsgin thesVamiable Length 


List Selection Process with Three Locations © 

In this section, we develop expressions for the 
probability of incurring a specified cost in response to an 
arbitrary demand. 

Le teatime Lb yebe thescosteo., themmth transfer under 
PicompGOlneveOtmoe lust sot length be fromecquations (3. 16)" i 


1s straightforward to show that 


mee) CL 
r- — b ; 


Pr{C(n;L)=0|/M(n-1)=(j,k)} = 1 = Db - ae 
ce Ca) 


an 





i 
= 
Es 
— 
= 
ak 
le 
Ww 


Pr {C(n;L)=C,,(M(n-1)=(j ,k)} = oF az b 
() (,) 
fonts) 
oJ M-3 k 
Pr{C(n;L)=C,,|M(n-1)=(j,k)} = ‘Le + ag - ee ; 
() (; ) (7) 
and 
M-k J jtk J 
C, ) (|) (aa) (; ) 
Pr{C(n3L)=Cac(M(n-1)=(j,kK)} = aR b- gb+t wig- Eg, 
= ()) C1) (1) 


Once equations (3.16) are solved, the above can be unconditioned 
and numerical results for the probabilities of incurring the 
various costs can be obtained. 

The above can be approximated in the steady-state case 
by letting f and h denote the steady-state fraction of red 


and black balls, respectively, in the urn M*. We then obtain 


Pr{C(™3L) = 0} = 1- (1-£)'r - (1-hytb - (£+h)g, 
Pri{C(™;L) = Cpp}= hr + £'b, 

(3.19) 
Pr{C(™;L) = Cpo}= fg + (1-f)"r - hr, 


and 


Gioye - £25 = ieaee B ee 


Eric (ose — Cagis 


The solution to equations (3.17), when substituted into the 
above expressions gives the approximation to the probabilities 
Pic CMincmamcos twine LeCSDOnSemtoa an arbitrary demand. 
Unfortunately we are unable to obtain analytical 
results to the above approximations. Thus in the three loca- 


tion problem we are unable to observe directly the change in 
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the problem with increases in L. The next section discusses 
some of the obstacles encountered when attempting to gener- 
alize to an arbitrary number of locations. 


6. The Ehrenfest Decision Model with More Than Three ~ 
Locations 


The preceding sections have shown that the extension 
of the selection process problem to more than two locations 
greatly increases the mathematical difficulties encountered. 
The resulting equations of the model rapidly become less and 
less tractable. 

One, of course, would desire to exhibit some rela- 
tionships for the cross-transfer rate c(») as a function of 
the length of the selection list (L), the number of records 
available to the personnel distribution center (M) and the 
number of locations in the AST rerio denoted N. Failing 
this, however, one might be satisfied with the intuition 
gained in the two location problem. 

For instance, in the two center problem c(@) was 


et) This 


Shown to be approximately proportional to 2 
rapid decay points to the increase of the selection list L 
as a means of reducing the costs of inter-location moves. 
The decision to increase L now seems intuitive in the 
MC te Ocat On case . 

Should analytical results be desired in the multi- 
location Casle, one first would desire to determine the scope 
OL Jules oulene We formulatesamnethod to describe the scope 


Tmo nomen in cnie Muiltir-Location Bhnenficst Decision 


Model. 





In the above sections, two methods of approach have 
predominated. Backward difference equations have been 
stated and solved. Also the state dependent Markovian 
properties have been exploited. General statements can be 
Pade alperne the scope of the analytical problem under both 


methods. 


Lemma 3.4. The number of conditional probability statements 


required to specify the system with N locations is N*-N+l1, 


Proof: There will be one equation giving the conditional 
probability of remaining in the same state. For each of 
the locations considered to be a source, there are N-1 lo- 
cations which are possible destinations. The conditional 
probability of a transfer from cone location to cach other 
is@the basiis)for one of the biackward equations. 

| Hence N(N-1) + 1 equations are required to specify 
the system. 

Note however, that the specification of the system 
equations is just a starting point under this method. The 
condition has to be removed from the conditional probability 
StacementSmand the result interpreted] se ingthe above sec- 
tions, we found the probability generating functions to be 
useful in the interpretation. it was through the probability 
Pole natliCmEMnctlonse that we were able to see thiat the com- 
position of urn M* obeys a binomial distribution in the two 


ieee ca anaed ‘Crinomidal’ in the three location case. 
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Now we consider the state dependent Markovian ap- 
proach to the description of the system. Let P(N;M) be the 
matrix of one-step transition probabilities. It is known 
that P(N;M) is a square matrix with the dimension being 
the Paneer of possible states in the system. Let |P(N;M)| 


denote the dimension of the matrix P(N;M). 


Lemma 3.5. The dimension of P(N;M), i.e. |P(N;M)|, is given 


M+N-1 
by Cy _y)- 


Proof. (By induction) 

It was shown in the two center problem that |P(2;M)| 
et tH so, eres) = Mess 1. “tienee the lemma ish true 
for N = 2. 

Suppose the lemma is true for N=n. We snow this 
implies the lemma is true for N=nt+l. First suppose that the 
n+ist location is excluded from consideration. Then the 
dimension of the one step transition matrix is |P(n;M)|. 
Suppose now that exactly one person from the nt+lst location 
1s in the urn M*. Then the remaining M-1 personnel in urn 
M* are distributed among the first n locations, thus giving 
rise to a one step transition matrix with dimension |P(n;M-1)|. 
aims) processwcontinues untzr!l thereg@are M personnel from the 
Te iste locatrom in the urn M* jand hence none trom any other 


location. Thus there are 
ferent) Mea (tah lo) | Sees (i; 0 ) | 


possible states in a system with ntl locations and urn sice 


M. Thus 


AAS, 





MRssc\s J. Male al: n-1 
ic ae st Ge 7 Jeo = (8) 


Observing from the Pascal -“Lriangle that 


ktn, _ k+tn-1 k+n-1 
CE) Ga i Geet 
k : 
: ane | 


it is readily seen that 


M+n 
[PGi isi) |RS* Ge ) 
: . F _ M+N-1 
Hence we conclude by induction that |P(N;M)]| = ( Nee) eis 


Naa 2 

Ticme | Cheiiomor tie Matrix fF (N 3M) ane palVvel oy ae 
most |P(N;M)|* probability seTISES Consider, however, 
Ghie nmatureyof the! allowed transitions. Due to the restric- 
tion that only two elements of the composition vector can 
change in response to any given demand, that the change 
incurred can be no larger than one unit, and that an in- 
crease in one element must be matched by a decrease in an- 
other element, considerably fewer probability statements 
need be determined than the maximum stated above. 

Once determined, the matrix P(N;M) raised to suc- 
cessive integral powers gives matrices approaching a matrix, 
the rows of which are the limiting steady-state probabilities. 
Thus if [P(N;M)| is not prohibitive, results can be deter- 
Nencemmeniso. == 1P can be solved for the limiting steady- 


state probabilities. 


0 





Thus we have been able to at least quantify the scope 
of the problem for the multi-location Ehrenfest Decision 
Model. The latter approach, i.e. using the state dependent 
Markovian properties, seems to be more conducive to the ob- 
taining of results than the probability generating function 
approach. Although the initial statement of the problem may 
be tedious, once stated the problem may be solved. 

Themicnptneorecne SCl6Ctilomelist entcersethe problem 
in the determination of the elements of P(N;M). Thus as L 
is varied, the probability statements forming the elements 
Cup (NeMe must be Changed = Since this 1s the portion of 
the problem which may be tedious, approximations may be use- 
mieeat thismestage. As in section C.4, it may prove helpful 
eer wee EecmpsODLem aS One With replacement.” iuis 
leads to statements that are in terms of the Lth power of 
Mmerractionwsracther than combinatorial®terms im‘. Thus the 
approximations may be more readily calculated than the exact 


expressions for varying L. 


D. DIFFUSION APPROXIMATIONS 
1. General Discussion 
Diffusion processes are Markov processes in which 
only continuous changes of states occur (sec Cox and Miller 
[7ljepaee 2035)5 For a more mathematical description, sce 
Feljer (| l¢ |MeeDittusition processcs are often@®derived by the 
scheme of dividing the time axis and the displacement axis 


(for saeOnce@incnsional process) into a very large number of 


IS) ll 





small intervals. A discrete Markov process, such as a ran- 
dom walk, is then observed to occur on the fine grid thus 
established. In many cases the limit of the discrete Markov 
process as the increments on the time axis and displacement 
axis aucerorced to become smaller and smaller as§a diffusion 
DILOCcCess!| =): 

gMore@itormally, consider aystochastic procesis 
Pete), athe Of. Let u (X(t) be a continuous, bounded func- 
iPmonmor X(t). SFurther, Met u(t,x) be tite expectation of 
ui (X(t)) Sonatetonal on the hypothesis that X(o0) Sex. milf 
u(t,xX) satisfies the diffusion equation. 


a, & Fook 
wit “ a 





where 


u(t,x) = Su, (ya, (x y)dy, 


q, (x,y) being) the density function for the transition 
pregaollity. eee tor the probabulity that the random tunc- 
PLO att iewilimoc in tthe interval (y.ytdy jmac, Gime ‘tf given 


Pihneriade ches valuemx at time Um and 
iPipmeboime1a xed A> Oy Liy-x|oadt O% YI dy 1S Oi ae 


= f (y-x)q, Cx, y) dy>b (x) ase tO, 


(2) 
ly-x|<A 


iirc 


ct] 


(3) (y-x)*q, Gy)dyra(x) as tod, 


ly-x| <A 


eo 





then u is said to be a diffusion process. The limiting quan- 
tity b(x) is often called the infinitesimal velocity or 
drift. The quantity a(x) is the infinitesimal variance. 

For more detail and examples, see Feller [10]. 

In other cases, a discrete Markov process may be ap- 
proximated by a diffusion process. This approach is essen- 
tially the reverse of the limiting scheme discussed above. 
The procedure is useful because of the resultant gain in 
mathematical tractability of the problem when the continuum 
1s used for the sample space. In many cases the use of 
diffusion approximations to discrete Markov processes has 
enabled the constructors of stochastic models to describe 
time-dependent phenomena that otherwise defied description 
Guest. obscrred the insi@ht provided by description 
through the use of more conventional methods. 

For instance, Gaver [12] applies the use of a dif- 
EUStOn fapproOximations to the description of the meniRtn time 
Wise he M/G/1 ductecemm ne themintroductiomeros that spaper eine 
States "the understanding of system performance furnished 
by the presient mathematical theory is) inadequate. The rea- 
son is that while the consideration of simple problems 
typically yields elegant mathematical results, the form of 
these results...is not immediately comprehensible nor 
useful in simple comparisons." 

Gaver and others have found diffusion approximations 


useful in many other cases [13,14]. Many other people have 
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used the diffusion approximation, including Newell [31,32], 
Kingman [25], Kimura [24] and Iglehart [22]. 

There are two basic approaches to the application 
of the diffusion approximation to a particular discrete 
Markov process. One, used by Gaver, Lehoczky, and Perka 
[15], goes directly to stochastic differential equations 
which describe the movements within the system. The solu- 
tion to that equation then leads to the approximate dis- 
tribution required. Another, given in Cox and Miller [7], 
involves an examination of the fundamentals of the limit- 
taking procedure. The limits must be taken in such a 
manner as to not obscure the randomness that is being 
J] vet liidt ts. Liste I1incrementsmonmrne displacement 
axis decrease as fast as or faster than the increments on 
enema cine axis, thesnresulting™ limits give awdetCerministic 
system. Also, as the increments on the displacement axis 
become smaller, the number of increments must increase in 
such a anor thatetne product Of the®number on imecrements 
and the size of the increments increases. 

Both of the above approaches have the same end 
result, so they appear equivalent. We choose the latter 
approach, since the limit-taking process provides insight 
Into. the application. 

eee liewlwo Location Problem with L=1 

This section devclops the diffusion approximation 

Pome nOdcimot section B.2 above. in ordér to proceed, it 


iS mec@ssary to rédcefine the systcm parameters. Although 
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in the previous development it was found to be sufficient 
to specify only the number of reds in the urn M*, here we 
use a function of both the number of red balls and the 
number of black balls in the urn. 
Consider the process Y(n)=R(n)-B(n), where R(n)+B(n)=M, 
the total number of balls in the urn M*. The process Y(n) 


MS characterized by 


pia + Lee (Tes (tele (3.203) 


Where Z(nt+l) is a random variable which is dependent on Y(n). 
Enes dustribution of Z(nt+l) ws developed tshortly. In keeping 
with the philosophy of the diffusion approximation, the 
magnitude of the step is to be variable. We assume there 
mea tGallsein the urn. Mf R({n) increases by 4, then B(n) 
decreases by A, and the change in Y(n) is seen to be 2A. 
cumilariv> 1i Rin) decreases*by A, then B(n) increases by A; 
and the change in Y(n) is -2A. We then have that R(n) + 
B(n) = MA. Then MAtY(n) = 2R(n) and MA-Y(n) = 2B(n). 

Adding MA to both sides of equation (3.20), using 


the above relations and dividing by 2, we obtain 


emt = 8 R (ne 257( lee 


ius wae characterization Of the process described in sec- 
MmiCimebeze te A—-l. Thus trom (53.)) we cam write dinectly 


EReeanEStributiton of Z(n+1): 





Pr{Z(n+1)=24|Y(n)} (es MA 


Pr{Z(n+1)=0 | Y(n)} 6: : 








DIR Ne 
+ 
Aes 
_ 
" 
m< 
Vee. 
3 
—— 
4 
! 
=: 


and 


|: 
ar 


Pr{Z(nt+1)=-2A|Y(n)}= ( 1 + 





a 

> 
= 

NO 


Now the expected value of Y(nt+tl1) is determined. Tak- 
ing the expected value of both sides of (3.20), conditioned 


on knowledge of Y(n), we obtain 


E{Y(n+1)|Y(n)} = A(2r-1) + (1 - @) Yq). 


Letting p=2r-1 and ael.$ we have E{Y(nt+1)]Y(n)} = Ap + aY(n). 


M? 
Given an initial condition Y(0) = R(0)-B(0) = Yo» we have 


EQY(1)) = u(1) = Ap + ay. 


eid O0 een Cire 1a Il 


E{Y(n)} = un) Ap (1tata2+...t+a™ fy + avy. (5229) 


Proof. The lemma was shown to hold for u(1) in the above 
Fr ecussi on The proof is then immediate by induction. 

Now for A=1, the limit as n grows large of E{Y(n)} = 
~P— = Mp = M(2r-1). But Y(n) = R(n) ~- B(n) = 2R(n) - M. 
Alsomt romagectionssB. 2, we have ER(~) = Mr. Hence, 
EY (n) = =e M(Z2r-) mee Thusgeche procctism, (nj—astat least 
consistent with the earlier work in expectation. 

Mime pacssion 1£0r themvariance of Yin), denoted 


VY eo r o“(n), 1s now required. We obtain 


Vecnemeiv(n)} = A-(1 + 4r(l-2)) ) = ¥(n)*/M. 
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When the condition on Y(n) is removed, we obtain 


o-tntljye=— (1 = =) eon) + 8A2r(1-r) 
(3725) 
Yo 2n ( Yo 


= a 2Ap ( = a Ap ) - a =e No ie 


Eomeeno ad nNiicrateconar tion X(Oj8—- yosewe have G-(0j)=] 0. bet 


O 

y 
Des a - c, c = 8A?r(1-r) and d = ~O - dp. Then (3.23) be- 
Comes 


Geant dee boo (in) + em a" 2Apd = ong 


For n=0, this becomes 
G7 (@)8=- c = 2Apd = d?. 


Lemma 3.7. In general, 


n Ca Zn 3 
1-b a -b 


G- (a) — ae. Cy) 








Proof. Again, the proof follows immediately by induction. 


Again, for A=1, o2(#) = ~S = eraErai = 4Mr(1-r). From 


section B.2, we obtained V{R(n)} = Mr(l1-r). Hence V{Y(n)} = 


V{2R(n)-M} = 4VarR(n) = 4Mr(1-r) is again consistent. 


Lemmas 3.6 and 3.7 give expressions in the discrete time, 
indexed by m- We need to change fromgthis time scale to 
apeowemmueus tame scale. We do thils by s¢ttiing t = nt and 
Petting, 1, the time interval between transactions, approach 


thy Oe 





For the diffusion approximation to be applicable, 
we need u(t), the continuous analog of u(n), to approach a 
constant. We also need the variance to approach a constant, 
i.e. o*(n) approaches some constant k. Note that if k were 
zero, the System would become deterministic. To avoid the 
possibility of the system becoming Gece 1 Mm Tis tae eaten TS 
necessary that the step size A decrease at a slower rate 
than t. The choice A = 17 1S appropriate. 


The relationship among the parameters of the problem 


are given: 


fe 3 yt’, 
As 14, 4s) 
MA? = k, 
and 
i ne. 


where y and k are constants. 
Substituting into equation (3.22), we obtain 


t/t 


u(t) = 2yk + (1 -—) 3 (y - 2yk) 


O 
which in the limit as t approaches zero becomes 

u(t) = 2yk + Oe = 2yK)exp (tk): (pS 2G) 
Some larity? Substatuting into (3.24) and allowing t to ap- 


Proagel zero, we obtain 


o*(t) = k ( 1 - exp(-2t/k) ). (32 2 
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In the limit as t becomes large, equations (3.26) 
and (3.27) assume the steady-state value. We obtain 
u = lim u(t) = 2yk 
tre 
and 


Se = mimeo Cale) kK. 
t+ 


Since k and y were arbitrary constants, the last relation 
determines k as the steady-state variance and thus y is seen 


to be u/202, Now (3.26) and (3.27) may be re-written as 


i ee Oy qm IEXp(7t/o5) oc) 
and 


He) 


o? (1 - exp(-2t/o7)). (5.299 


Having obtained expressions for u(t) andd(t), we 
now develop the distribution function which Y(t) obeys. We 
can state that 


(teem (te. 2 (Cit) 5 


where the distribution of Z(t+t) 1s obtained from equations 
(3.21) by direct substitution. Looking then at the moment 


McheTatiiiy, Hunetron Condit LOnca on, Y(t Je 


icra ee — Etexp (mY (irr lmpe be edexp(m:! (ite) )|Y(t) }. 


% c mt? mt? 
mt -MT = : 
Letting A=e = and B =e +e , we obtain 


ie 
O(tstscm) = O{tem) + 4O(t ;m)ACEZytT “B) 


(t : ee 
- : 1 oo Ct sm) A(B+2yt 7A) ) 
4Mt 2 1. 
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which in the limit as t approaches zero becomes 


20(em) at « zymocemy ~ cB 20ftim) 


We now have a partial differential equation in two indepen- 
dent variables. Substituting y = u/2o0? into the above equa- 
tion gives 


a onee, m1) 


{(u + mo*)@(t;m) = 


o0(tim 7 em 
oe = 1 ie are o108) 


o 2 


To solve this partial differential equation, we follow 
the procedure used by Cox and Miller [7]. Let K(t;m) be the 


cumulant generating function, i.e. 


"Ket -m)e2e an 6 (tm). 


Then equation (3.29) becomes 


9K (t 3m) , _m aK(t3m) 


“y) - 


With the initial condition Y(0) = Yo? equation (3.30) has 


the solution 


K(t3;m) = mu + m(y -H)exp (-t/o*) + Le (l-exp(-2t/o7)), 





Cmmeoe ec omenasmence solution exp(K(t- mig eenomeyhich (Ft ;m) 
ms @readily yecognized as the*™generating*™ function for the 
Rom cimeals erp outIon. FurthermoregsUusimemeqmations (3.25) and 


(3.26), we have 


o(tim) = exp(mi(t) + me BG). (3.32) 





Hence Y(t) is seen to be a time varying normal process with 
mean u(t) and variance ¢(t). Further, as time increases the 
process settles down to a stationary normal process with 
panameters w! and o>. That ils, 

lim @(t;m) = 9(@3m) = exp(mu + 4m*o7). 

to 

Thus we have shown that the selection process for 
the two location institution can be approximated by a dif- 
fusion process. From this approximation, one is able to 
obtain insight into the nature of the transient process, 
1.€. how the composition of urn M* behaves following some 
initial disturbance away from the steady state case. 

In the next section, we develop the diffusion ap- 
proximation for the two location, variable length selection 
fust problem. 


3. The Two Location Problem with Variable Length 


selection List 

Die tiie seet1on, We Extendethesresults or the pre- 
vious section to the case of the variable length selection 
list. Again, we examine the process Y(n) = R(n) - B(n), 
HNeCuCMeNCEprocess i5@eCcharacterized byeequation (3.20). 

In order to define the conditsjonal distribution of 
Cine Deveney (Mjieett 1S)convennient £6 wsieythe approximation 
Ob section Be4. Using Equation (3.3a) from that section, 


the analog of equations (3.21) becomes 
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pr{z(n+1)=24|Y(n)} = 24a - SL) (3. 33) 
L L 
pr{z(n+1)=0 |¥(n)} = 1 - 27¥¢ca - SYD) re caeE My (1-1) }, 
and — 
L 
Pr{Z(n+1)=-2A|Y(n)} = Gd + ain). el rie 


Making=use of the rellations’ (3.25), the development 
is quite similar to that used in the previous section. Let 
MenuneaciOtcmtne cexpceted valuc Of 1(t) Under the policy 
Ofmumscuection dist of jeneth L. Aliso let o*(t;L) denote 
Ghewvarianee of Y(t) under the policy of a lisit or Wength 
Faeroe olvelicnt an this €ase €o amkefthe arbitrary 


Polstait y Bemastunction of b, dénotedsey(L) where 








yibe= Y(l)ie y(1) being the y of the previous section. then 
his 
u(tsL) =u + (yj-wexp(---= ) (3.34) 
2 ei 
and 
; 4Lt 
oo (tbe = —- (1 > Geen = = ; ices ige (5.5515) 
Oo 


where u and o? are the limiting results obtained in the 
mueececeuimee SCetION. 

The generating function is now also indcxed on the 
HeTmc Otetnecm list &, and 15 denotcd Gfe .m,L).) It as rela- 
tively straightforward to show that 


2 2 
Ge ney) e exp Cun Gagan (ttle, 





Witgwens again recopnized to be that of a normal process. 


eminem icttieac eh increases, the variance of The approximation 
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goes down by iL as steady state is approached. Also, the 
exponential decay rate decreases by a factor of We 
Thus, there will be less variance once steady state is 
achieved, but it will take longer to attain steady-state 
conditions. This 1s intuitive, since the selection process 
described by the variable L case is tending to maintain the 
Status quo, i.e., is designed to reduce the flows between 


locations. This in turn prevents the system from rapidly 


attaining steady state. 


Ege LENSIONS AND CONCLUSIONS 

bieulimmendapter sswe Mave developed wa model ofethe person 
nel selection process within a large institution. The selec- 
tion process was seen to be first concerned with finding a 
Suitable (qualified) individual in response to demand for a 
penson with Specified gqualiffications. Of secondary concern 
Nase tie Current location of the) selectees, 

As the length of the selection list, L, was allowed to 
grow, we saw the steady-state inter-location transfer cost 
decrease. Thus, if the institution 1S concerned about de- 
creasing the transfer costs then the personnel managers 
should be tasked with finding more qualified personnel prior 
Eome ut cetin ved Pr ans terse Unde Pere meccUNDEMenom os the chapter 

The model was developed for an institution with two lo- 
Cations, then extended to the case of an institution with 
thmeemilocations. Comments about extensions to institutions 


With more locations were madc. 
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Approximations were developed for steady state results. 
Also diffusion approximations were developed to aid in 
describing the transient flow given some deviation from 
Steady-state as an initial condition. 

Throughout the chapter, we were implicitly concerned 
with the minimization of the transfer cost in response to 
individual demands. Such a myopic view of the selection 
process does not consider the long range needs of the in- 
stitution, nor does it consider that moves effected early 
in the period can influence llater moves. That is, of the 
L selectees the one associated with the least transfer 
G€osii may besmore Suitably held in thesurn M* for @ later 
demand. In that case, it might be desired to transfer an 
individual with a slightly higher cost of relocation. 

ie nosilts  etechapter fl arepsecn to mimpinge scmermec 
personnel selection process at this point. The use of the 
model as a personnel planning model gave as an output the 
Sptwmal sleCworsinter-locaticm) transters by pay grades If 
pay grade is the determining factor in the qualification 
Specmpicdsimetamspechapter, then the amstitution would de- 
Siremeomcimcct a wtranster only af thats particular trans#er 
were in the optimal set provided by the periodic planning 
model. Of course, this may not be possible, since it is 
hommelanantcad that gucheagtransfer 1s inelmded in the 
options upon completion of a particular selection of L can- 
didauesmeouch an event would imply @ higher cosit of tranis- 


fers than was provided by the planning model. lIlowever, the 
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use of such a scheme would tend to avoid the gross subopti- 
mization of the selection process without such considerations. 

In the next chapter, we investigate the problems encoun- 
tered by the personnel distribution center when it is forced 
to select more and more candidates Dimon tLOnetiectrmg aa 


transfer, reflecting an increase in L. 





IV. THE RECORD SEARCH PROBLEM 


A. GENERAL DISCUSSION 

In chapter III, we saw how the expected cost of transfers 
1S reduced by increasing the length of the selection list. 
In this chapter, we develop a method of determining the 
Lime wreguired for thle personnel distribution center, M*, to 
establish a selection list of length L. Obviously, if the 
time required to establish each individual list is too great, 
the personnel distribution center could encounter operating 
duiticulties. 

Let us approach the problem from the standpoint of the 
Pelsonne! distribution center. Suppose there)yare’ a ehixed 
number of records, representing personnel available for 
tEranster, wit theldisposial of the distribution center, Let 
the number of records be denoted N. These records may be 
considered to be the personnel files for individuals physi- 
cally located at many different IGEMEONS. From the N 
records, the distribution center M* must either 

(1) select an individual to fill a vacancy (L=1), 

(2 iy se Vee ema woroup if individuals as candidates for 
tne Vacaneymivarlable L, Ligreater than 1) or 

(3) demonstrate the there vs no qualified replacement 
among "the? N available records. 

The procedure to be followed will be an investigation 
Cie ticmproebLven wien the number of records 1s véry large. 


Restrictions placed on the system iby increasing the length 
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of the list L and by restricting the number of records 


available are then investigated. 


B. THE BASIC PROBLEM -- UNRESTRICTED NUMBER OF RECORDS 
1. General Distributional Results 

Suppose that each record made available to the dis- 
Eripvlcnonmecenter M* Was a sprobability m(@ji) of beimeythat of 
a person qualified for a job of category j. For brevity, 
we will say that this person is of category j. Suppose 
m(j) ts’ independent of all factors*except the category 
esiel i 

Let M(j) be the number of records among the N avail- 
able which represent personnel of category j. Under the 
independence assumption above, M(j) is distributed binomially 


with parameters N and m(j), l.e. 


Pr{m(j)=n} = (4) m(j)"(1-m(j))*™. 


PEevoe cine the number jon records Otmpersonmernoer 
of category j searched prior to finding the kth person of 
category j. We assume that the records are observed in a 
random fashion, i.e. that no sorting has been done or any 
other procedure which would tend to make the qualifications 
of one record dependent on the records searched in the past. 
We ignore momentarily any restriction on the number of 
records available. That is, Wi N Ws large enough, we will 
always) be able to find the k personnel having the specified 
(HUjtpionedt rane i nem, '();K) 1s#seen to have a negative 


Pons tmiiii1 Oo With parameters K-and m(j). That is, 
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Pen het C mG) C-m(j) 0". 


Suppose now that we restrict ourselves to searching for 
personnel of category j}. We then can suppress the j, 
letting m(j) =m. 

Given F(j;k), the total number of records searched 
CO Liles che set person of Category 7 is then F(j;kK) +%#k = _G(j 5k). 


Let us define the generating function of G(j;k) as 


C(j kez =] Biz? >}, 


Then 

~ Te ° = Z™m 

6(53k32) =(q57y) - (4.1) 
Now the expected value of G(j;k) is given by the first 


derivative ot (4.1) with respect to z, the result evaluated 


at z=1. This gives us 
E{G(j;k)} = k/m 


Also, the variance of G(j;k) is given by 


mn 


We a 


Saige + REMC mee (Gi(jgk)), 
1 


Z= 


from which 


k(1=-m) 


V{iG(jsk)} = 5 


We thus have obtained an expression for the generat- 
TeoetUNCtLONmoE tlie number of records Searched and have de- 


PeTMined= tie mean ana Variance ot the distribution. 
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Again, leaving N, the number of available records, 
unrestricted, we introduce the second attribute of the 
selected candidates. Here the second attribute of concern 
is current physical location. Let &£(j3;1) be the probability 
that an arbitrarily selected person of category j is cur- 
rently@located at 1. Suppose &(33;1) *to.be¥a function of 
Salyecme category af and’the location ti.) Let 0Of7:;1) ibe the 
number of qualified personnel of category j selected in 


Srdcimremontaim the first person of catepory 7) located at wa. 


I 


Pr{Q(j,i)=k} = (1-2(j,1))°12Gj,i), ke=1,... 


Also 
OC) 51) 1 = Z&(j 51) 
Biz Ss eee oe 
Sr jemmeisescen ago be distributed eecometricallly with 
parameter %&(3,1,). 


Nowslet R(},1) be the number of records ysearched to 


fiideche first qualified person of category j at location 1. 


Then | 

BE { 28Ci-2) 1a¢j,ijp=k} = BE { 260 >K)}, 
Witeam1Secaven by (4.1). Thus, removings the condition on 
Q(j,1), 


Cas) ee zm(j)%(j,1 
E8009) = STG DT (4.2) 


FROMewAtenenN|jed) 1S)'Sé€cn to be distributed geometrically 


Wheweparameer m(j)2(j},1). 
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In the next section, we discuss the distribution of 
the time required to search through the R(j,i) records to 
find the first qualified person at”’the specified location. 

Ce LINCmRCUlarements tor Basic Record Search Problem 

Given the number of records required to be searched 
to find the first person of category j at location i, we are 
interested in determining the distribution of the time 
expended in the search. 

Suppose the time to search the kth record is denoted 
X(k,t). Suppose all searches of records obey the same dis- 
tribution and that furthermore, the time to search any 
pecolGmieomimdcnendcnt Ob the time to Searcheall other re¢e= 
ords. That is, the times are independent and identically 
distributed (i1.i1.d.) random variables. 

Lop] i bcm the time =leduiirca= to Lind eme temas L 


Per sOMmOLTerkctOoGmumy dt LOCdtion 1. tien, for dee iyen K(y i je 
Tels ls) MO EOE men esr GONG PS ae hey 


Sineewn( j,2)eus the number of records which must be searched. 
binge pthes random) variables ares if isd. @thcek#can bemsuppressed, 


and the random variable is now denoted X(t). Now let 


X(s) = E { SPC 

Then 
E { a! Os) lcqeeion } = X(s)", 

from which, by removing the condition on R(j,i) we obtain 
T(j,iss) = Bf @ SPU) 
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P wie 
a aces (COLECLIGHETG ES OPA (eo 8) 


which is easily shown by substitution into (4.2). From (4.3), 


‘ 


it is obtained that 
ET(j,i) = ER(j,i)EX(t) = EX(t) 
and 
VarT(j,i) = ER(j,i)VarX(t) + B*X(t)VarR(j,i) 


~ Var x (en 2 1-m2 
ae 7a Eo. Ge Tmty2° 


In the special case that X(t) is the exponential 
distribution with mean nee, T(j,i1) can be shown to obey the 
exponential distribution with mean fum(j)2(j,i)} 72. 

Thus under assumption of independent, identically 
distributed times to search individual records, of an un- 
limited number of records available, and of no restrictions 
SMmrnce number ot Qualified personnel @found prior to effecting 
a transfer, we have an expression for the transform of the 
@ustribugion Of the time to find the finst peyson of cate- 
gory j who is currently located at 1. If these assumptions 
heldnmencmsy stemmecould be operated wilh mo Gross-transfers 
and hence no inter-location moving costs, provided only that 
tnemeime to@S@ert through the possibly large number of records 
1s immaterial. 

We next look more closely at the reason for search- 


ane to find w number of qualified personnel prior to making 





a decision as to which person to transfer in response to a 


given demand. 
x 


eRe otelCnlots Ove ine LENGHT OF THESSELECTION LIST 

The requirement for the personnel distribution center 
to search until finding L qualified candidates as was re- 
quired an chapter HIT, ws) somewhat arbitrary. In thils sec- 
tion, we investigate an alternative formulation of the 
problem. 

Consider the two location problem of section III.B.3 
where the variable length selection list problem was intro- 
duced. There, after selecting L candidates in response to 
a demand, the current locations of the candidates were ex- 
Hi LLeCmeomUetermiane Whether Gr Not any one of the candidates 
were at the location from which the demand originated. 
Another policy might be to examine the current location of 
eachmse Mected) Candidate to determine it he is¥at the loca- 
ClonmEesOn Whremetne demand oOrnlginatiea.)) Ine Search continues 
Until either fhe firsit person at the proper location is 
selected or L qualified personnel have’ been selected, none 
of whom are currently at the location from which the demand 
Srroinatcammmmae this CaS€s m5 1S aAheuUppecmelim?eeon thes number 
Smgudlitica repliacemenus to be fond prior, to effectingla 
Cldil> [Cm smi Clhesprevious case, I) istamcontrol variables 
Wwe. Subject to institutional policy. 

Suppose now that the length of the selection list is 


Bectpictcusto dmmaximum valuc, denoted L. Let R(j,a;L) be 


nay 





the number of records searched to find the first person of 
category j at location 1 or L people of category j, which- 
ever occurs first. 

Since the attributes of category (or qualification) and 
location of the individual represented by any particular 
record are assumed independent of any previous records 


searched, a renewal theory argument tells us 


Rigieesl) = el with probability m2, 
= 1 + R(j,i;L-1) with probability m(1-2), 
and eR eds) oWwethmrobabiinity l=me 


where m = m(j) is the probability that an arbitrarily selected 
He COlrd 1S wommeatcmorven, £ =f (1351) ws the probability that 
an arbitrarily selected person of category j is at location 
Pere eee ts 2 random variable with tht@same distri- 
bution als, R(j,1;L). The relations hold for L greater than 


one. For L=1, we have 


Raa) il with probability m 


1 + R'(j,i3;1) with probability l-n, 


and 
Simeceror L-leeehnessearch terminates’ upon selection of the 
fie Sid Ua Ll ieveda person . 
Hence sities cenenating funety on R(z:1) = ih pe oT 1s 
given by 
fled) = eases fe (4.4) 


Wencal alse casily obtain theseecnerating function 


RQ(z3;L) = B(z8G 21:55), eC VismeOnen ezine line ercSsult as 


-153 








R(z;L) = g 


Hii 3 (1-2)R( rid a 
a Z (a5) 


GR 
Defining R(z;0) = 1, the above reduces to (4.4). 


Lemma 4.1. The generating function for the process with a 


selection list of length L is given by 


mS L 
i 1-2 
a Toad) (2 + (1-2)(7 tiny ) ).4.0 
PrOG@mn (MeN prooL 1S immediate by induction. 


From (4.6), we obtain the limiting result, 


R(z3@) = lim eG Ze1e3) Sa TaD 9 (4.7) 


L+< 


which is the result (4.2), obtained for an unrestricted L. 

From (4.4), we observe that for L=1, the number of 
records searched is geometrically distributed with parameter 
m. For unrestricted L, the relation (4.7) is observed to be 
Les tierra tINQ@erunctlOn Ofederandom= Vartabilc Obevyimy the 
geometric distribution with parameter m&. Thus as & in- 
Srecascs melee sty 1bDUtLION of Rij ,1 Emo ocsetrom peometric in 
[muomrcometricmin me. 

ite _ of interest to observe how the expectation 
Ome GeaMevaries with lL, @@this mayeebe achieved 
DyeESUleCCISi Ves dlitverentiatlon of (4.0). 

By differentiating (4.6) with respect to z and evaluating 


tite result temiaz-1, we obtain 


L 
_ ween 
E{R(j,i:L)} = 2 aad (4.8) 
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Equation (4.8) can be used to give rone insight into how one 
might choose the length of the list. Let R*(j,i) be the 
number of records that the personnel managers "expect" to 
search in finding a person of category j at location i. 
Suppose also that m and & are known for all locations and 
categories. Then, equating R*(j,i) to the right-hand-side 


of (4.8) we have 


iF 
‘ae 1 - (1-2 
met) 3S Se 


Or 


p= Mn CemeR*C5 4) (4.9) 


Noweeieek (7,1) aS on thle order of (me) 1, then we see that 


foese Vvetyelargesss in this, Case the’ results of section B 
above are applicable. 

On the other hand, if R*(j,i) is on the order of mt, 
then L is close to one, the minimum that L can take on in 
she Gentextwes the problem. “WE R*(j i is®lless! than Ti 
Ehen one would not “expect” to find any qualified personnel 
of category j in a scarch of R*(j,1) records, much less one 
who "rs Gurrently at a'’specific location: @hlence 
mt ek (jae) < (m2) is the range of intcrest. 


The next section is a discussion of the problem under 


hBestrictions oneN, the numbcr of records available to M®. 


f 


D. RESTRICTIONS ON THE NUMBER OF RECORDS AVAILABLE 

Suppose) mow there are a limited number of records avail- 
doLémuouthnempersonnce! distribution center M* from which the 
Samecetlons are te be inade. Then the introductory remarks in 
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section B give M(j), the number of records of category j, 
as being a binomial random variable with distribution 
parameter m(j). 

Also equation (4.1) gives the generating function for 
G(j;k), the number of records searched to find the kth per- 
son of category j. However, (4.1) was derivéd assuming 
that the number of records available, N, is unrestricted. 
To determine the effect of restricting N, let h(j;k|N) be 
the probability that there are at least k persons of cate- 


gory j in the N available records. That is, 


h(j;k|N) = Pr{G(j;k) < N } 
N . 
= nek Pr{G(j;k) = n}. (4.10) 
Kare 
= 1- nz0 Pir exactly me personsmor 
Category jain the N 
records } 
leail 
- 7 Dre O ac gee afl del 
= ] nz0 nd oie = M1} 


Cc Cem be Lic eNUmDemmOlmnCCOn dsm sCdrcled an. rec - 
sponse to a demand from location i for a person of category 
[ecu vOscmmiidemene "Scarech procecdssunti ls the first person 
Of Gategory } who 1s) currently available at location 1 ‘1's 
PoOuntmonrsaliaNerccords hayegbeen scarched. If all N rec- 
Oudsmare: Searched without findangy a person of category j, 
the demand ils £11@d into a’delay system to be serviced at a 


much later timc, thus ensuring independent searches. If 


Here 





Upon searching the N records, one or more personnel of cate- 
gory j have found, but none are at location i, one of 


those found is transferred, ending the search. 


Then we have 


R(j,13N) 


tt 


G(j,k) ue the kth person found 
is the first to be located 
at i and G(j,k) < N, 


= N mi at) licast one person of 
category j was found, and 
none of those found were 
located it i, 


NH R'(,13N) ii ones Personuneilo. —Catero ry 
] Were found sn the sscanch 
Si wieeecords . 


where R'(j,1;N) is a random variable with the same distribu- 
tion as R(j,1i;N). Then the generating function can be shown 


to be 
Pej iiNee Bige de hy 


a © Soe IE ICED y 
hoe e Tz (1 -m) Ih= A(C sin) 


24 


2! ites) - am) 


It may be that the restriction on N is immaterial. 
Recall that R(j,i) is the number of records searched to 
find the figse person of cat®gory j who wsecurrently at 
T@oGation 1, where the selection list is unrestricted in 


Fenoclieeerol (deel c 15 Stralghttorward to show that 


ee 





il 


E{R(j,i}} (m2)- 


and 


VIR(j,i)} = (1-m2) (m2) 72. 


Now if N is greater than E{R(j,i)} plus some number of 
Standard deviations, one would "usually" expect to find a 
qualified replacement for the demand at the location from 


which the demand originated. That is, if 


it 


N > (mz) + y(1-m2) 72)” (4.11) 


for all categories j and locations i, then the restrictions 
on the number of records available to M* is immaterial. 


For instance, if R(j,1) is approximately normal, then 








Pr#R(j ,1) 


iV 


hae Pri Rk (jee 


ia Gd aia = ee 


and the latter term has values given by 





For a given N satisfying (4.11) for a given y, we see that 
the probability of not being able to locate an individual 
of a fiven category at a given location decreases rapidly 
Wateoeeatiimemedse an y. That is, the movements between lo- 


Greroiterare mot. Usiialliy’ caused by the restrictions on N. 
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Note that y is a decision variable. Given m(j) and 2(j,i) 
for all categories j and locations i, one can find N as a 
function of y by treating equation (4.11) as an equality. 
Doing this for all category/location combinations provides 
us with a set of N's. The element of this set with the 
greatest magnitude represents the number of records which 
must be made available to M* in order to ensure that the 
personnel distribution center is (to the degree represented 
by the magnitude of y) unrestricted by N. Increasing y im- 
plies increasing N*, the maximum of all the N generated. 
This reflects a greater requirement for M* to be able to 
lGcalcmasreplacement of any Category at any location. 

The multiplier y may be expressed as a function of loca- 
tion i and category j for all i and j. The multiplier wouid 
be larger for those locations and categories for which moves 
into or trom the positions in question have high costs rela- 
tive to other positions. For instance, in the three loca- 
tion problem of chapter III, we had that moves into and 
from location R were more expensive than other types of 
NG Comms Lie. sease. we might ssct thesmuleaplucr y to bema 
higher value than in other cases. 

Mahitemthne multiplaer ¥ agfunction of Wecation 1 and 
Gaeetom esp enassother applications. ~it may turn out that N*, 
Cre smaximumeet thegciatlculated N, 15 too larsesgfor reality. 
They valucseN(j,1) which make up the set of N firom which N* 
1s the maximum may have only a few values which are "unrea- 


cole inemy stor those particular Jocations/grades 
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might then be decreased, reflecting the decision to absorb 
the costs ery to be incurred due to the inability to 
find that particular combination of category and location 
in a reduced number of records available, and thereby re- 
ducing the value of N*, 

The determination of N* may govern the size of the 
personnel distribution center M*. If N* is large, then 
many personnel may be required to maintain the selection 
process. 

The next section formulates the problem for a restricted 
number of records available, N, and a restricted length 
scleeraon | iste el, 

Pee JOLNIe RESTRICTIONS ON THE LENGTH OF THE SELECTION LIST 

AND THE NUMBER OF RECORDS AVAILABLE 

The most general case of the record search problem 
includes restrictions on both the number of records avail- 
able to the personnel distribution center and the length of 
selection list that the personnel distribution center is 
mould. Gomestabmishepriomevo effectinggay@nansfer- 

It was noted previously that L, the length of the list, 
1s a control variable, fe. - subject to institutional policy. 
The number of records available to the distribution center, 
denoted N, may also * a control Vania Im ore instance, 


tic iruc Oienmeould be set by allowing only tiresnext N 


peoputeesciedulcd to rotate to be available] However, thiis 
Ways De edbtittcial, A*more realistic situation is now 
dee ra ied: 
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Suppose that all tour lengths for personnel of category 
j, denoted ie are fixed. Suppose also that Ne denotes the 
number of personnel of category j in the system. Suppose 
now that the institution sets an "availability period" 
which gives the period just prior to the scheduled rotation 
date during which the individual concerned is available to 
the personnel distribution center. Denote the availability 
period for a person of category j as ae Then the number 
Of peoplewor category j available to the distribution center, 


denoted Hho: 1s, on the “ayerage" 


a. 
ae = LN 
J - J 


Thus, the number of availabilities is seen to be a function 
of two policy decisions, the setting of tour lengths and the 
setting of the availability periods. 

Prloretosembanrking,on thesmathematicalmdescription of 
the model under joint restrictions on N and L, we discuss 
the sequential nature of the record search problem. Essen- 
tially, the personnel distribution center tsearches the@avalil- 
able records in response to a demand until one of four 
possible situations 15 encountered: 

(1) During the isearch, a person of category j at loca- 
tion isi found. The first person so discovered is transfer- 
Ted@andsthe@cost of searching thle*® records 1S incurred. 

(2) eine Search terminates with |, people in catepory j, 
NeVCeormniOlmicecirrently at location i, One of the qualified 
De sSonveimusmtnainsrerrea, and both a transf@r cost and a cost 
of searching is incurred. 
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(3) The search terminates after searching N records, at 
least one of whom is of category j. Again, a transfer cost 
is incurred as well as a record search cost. Note that the 
number of selected individuals in this case is strictly less 
than L, there will be less flexibility in the decision of 
whom to transfer, so one might expect the transfer cost to 
be greater. 

(4) The search terminates after searching N records, 
none of whom ois of category }. The cost of searching N 
records is incurred. Also, a "penalty" cost of not being 
able to locate a suitable replacement in response to a de- 
mand may be incurred. In any case, some time later the 
distribution center will again search for such a replacement, 
incurrink jaddiittionai costs. 

The sequential nature of the problem under joint restric- 
tions OneN ance is illustrated in figure’4. 1. 

Suppose now that the personnel distribution center files 
an unfilled demand into a delay system, from which the demand 
Wit) emMercemdamanlater time. Suppose the’dellay 1s so great 
Liaieeence personnel avallable ta the distribution center 
diUmainey Cle first search are Completely diritcremmumeron: those 
available during the second, or subsequent, search. Thus, 
the searches are micnendont: Pete Gee lenmenc =e tie enumber 
Of records searched to satilsfy a demand from location i for 
aeperson Ofcatecory 7 und@®r restrictions om N and L. Then 

i@edesneuem=  G{)] 48) if the kth person of catcgory j is 


Cie tio GGOmucmlOcated at 1, 
heer ee cticewani ges) <N, 
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Yes 













Transfer; 


Is indivi 
s individual Note 2 


at Gea ciom 
1 e 


Yes 
Iselections 
been made? 





Is? record of 
Gare CO ry 4 4 


__ Demand j/ Examine 
Set k=1 Record k 





No Yes 
Set k=k+l 











| Has at least one person | 
of category j been found? 


Noemie ineimecost of necordgsearch. 


Wo WCw7e, INGUr COSt of record#search 
Diliseeemansier cost; 


NOtcmereincir cost of recordssearch 
plus penalty cost. 


- Figure 4.1. The Record Search Problem with Restrictions on 
N and L, Showing Response to a Demand From 
hocation a for a VPersonegol Category j. 
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RG EN SL) = G(j;5) if the Lth person of 
: category j to be found 
1s not at location i 
and G(j;L) < N, 


= N if L persons of category 
j were not found, at 
least one qualified person 
was found, and none of 
those found were at i, 


Niece eg. 3 N 4) if no individual of cate- 
gory j was found in a 
search of the N records. 


The generating function for R(j,i;N,L) can now be found. 


We obtain 


RCP LON) = Teele Saat es 


i ll 


\f . 
= 2, E{z? 938) yn¢j sk [NJ (1-2) 9772 


me a a ee 


+ 2X th(j,1|N)-h(j,L|N) }Pr {none of those found 
were at location i} 


EE ich eto) 2) pie tical Noe 


After some rather messy algebra, we obtain 


zm& 


2a N N 
Gee eee Lez lem) i= T-z(i-mhy 


+ eet.) ( Aa ) 


N-n 


|) oI 
igs 


CN) m"(4-m) 


zm L- 1 


i Gas ) tema 
Tera, Lace (i-z(iem) 
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a eel oe (S) m(1-2) Ny-myN-2, 


=. to 
tt ot 
om) a 


Thus the general problem of restricted N and L does not seem 
to possess sufficient tractability to proceed analytically. 
ieee (acc mmOmeLneCmreastttetions On Other aspects Of thegdis- 
tribution problem remains a problem for future reconcilia- 
tion. Although there is not much insight to be gained from 
the mathematics, one does gain appreciation of the difficul- 
ties to be encountered in the analysis of the operation of 


thes personneledistribution center. 


F, EXTENSIONS AND CONCLUSIONS 

Mmciomenapter, we found thatwincreasing the length of 
the selection list and/or restricting the number of records 
available to the personnel distribution center affects the 
numbers of records that must be searched prior to making a 
deéeision te aiitectmastransfer. This in turn affects the 
amount Of time’ spent in searching in response to each in- 
dividual demand. 

MotiMemehem results of Chapter Tleimpht affectgthe search 
DeOCedUie mam OT instance, is taidmoursecareminos tor the first 
qualified person to be found who is currently at the location 
from which the demand originated, one might only search un- 
Br inciiemenes tirst qualified persom whose transfer to fill 
a demand is among those specified as being among the optimal 


Seo eethananers Under the provisions of chapter II. 
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Alternatively, one might set bounds on the number of 
Candi.aaces! coe perselected prior to effecting ja transfer. 
That is, there may be a value L denoting the number of 
CanatuatreS sOube selected prior) to erfecting al transter, 
provided no candidate is a "no-cost" transfer, i.e. is at 
the location from which the demand originated. Once L 
candidates have been selected, then the candidates are ex- 
amined to determine if there are any whose transfer would be 
in the optimal set provided by chapter II. If so, then one 
of those is transferred. If not, the search continues until 
either a candidate whose transfer would be in the optimal 
set is selected or L, a maximum number of candidates, is 
achieved. In the later case, a “non-optimal” transter is 
effected, The analysis of such a procedure would be much 
more complex than that presented herein. Since we found the 
EeciMMiEomOtoiiomenapter Co be difficult tovintenpret., and 
such extensions would be more so, we did not pursue the 
extensions. 

Another extension, also not pursued for the same reason, 
VOene MenCmCKanInation OL gehe CrLtectson allowing degrees 
Seequalitvedwc tom cOmecnucr the proplemsein Chis extension; 
personnel might be classified as qualified, over-qualified, 
marginally-qualified, and not qualified. The inability to 
Serccteaieiitim raul that towqualit#ed might result in the 
assignment of an over-or marginally-qualified individual. 
Such a problem would be faced with the determination of 


Mpenalty™ costs for using an over-or marginally-qualified 
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person in a position as well as the transfer costs associated. 
Here again, the concepts of a lower and upper bound on the 


length of the selection list migiit be useful. 
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V. CONCLUSION 


In this thesis, we have developed three analytical models 
which serve to describe rotations in multi-location, multi- 
grade personnel systems. The models, although basically 
independent, have been shown to interact on a macroscopic 
ey eqs 

The basic balance equations which we developed in chapter 
II are the first known analytical relationships describing 
thewcomp tex Interactions between ballet structure, rotation 
structure and promotion structure in a large institution, 
Under the assumptions of chapter II, the structures are seen 
to be heavily dependent on each other. Thus we were able to 
Show that the specification of the billet structure and the 
rotation structure completely determined the promotion struc- 
CU nee ii Smilmotndeorce of dependencegmay explain the probe 
Tencmencountercamim a personnel systempwhere the three 
STIrUCCUTCSmamemspecitied iae pendeunine 

In the development of the sensitivity analysis on the 
transportation problem, we found that changes in one element 
of the basic structures had complex ramifications on the 
constraint vector of the basic transportation problem. Thus 
we developed the "secondary" analysis given in section [II.I, 
where the feasibility of the primal problem was considered 
explicitly in the computation of the shadow prices. 

In chapter ill, we developed a model of the personnel 


Seer role yt @Gecsewithin a Large institution, where the 
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selections are made to determine which individual transfer 
iseto be madesin response to a demand for a’ person of a 
specified qualification. The model, which we called the 
Ehrenfest Decision Model, was developed for institutions 
with two and three locations. The model serves to parame- 
batizeaune spermsonnel=seliection process. The results spe- 
cify that lengthening the selection list, i.e. increasing 
the number of qualified personnel to be found prior to 
making a decision, serves to decrease the numbers of inter- 
location transfers which occur. 

In chapter IV, we developed a model which describes the 
problems encountered by the personnel manager when the se- 
lection list is lengthened. The results of this model in- 
PrcacCMmtlilatwones would mot desire to force the selection list 
to be arbitrarily long, since such action implies that the 
number of records Wnichmone ECxpeetSuuomme SCaTChedmm Order 
to specity the list may become very large. The large number 
of records would in turn imply a long search time and thus 
an increasing degree of unresponsiveness to arriving demands. 

im tChemeonclustons to both chapters ieand IV, the in- 
rlucmce SO weiemresulis of chaptcreWi yom the; models of those 
chapters was discussed. 

Finally, the pursuit of better and better understanding 
of the complex interactions between components of the per- 
Somme! systems of large institutions has becn our goal. The 
complexity 1s of increasing concern as more and more, institu- 


Piotemenretmenemtamily of mSs%etiitions with multi-location, 


169 





multi-grade personnel systems. Also, those institutions 
already in the "family" are interested in becoming more 

ana mone etticient in their utilization of personnel. Dis- 
economies caused by failure to understand these interactions 
are of concern to both institutional directors and their 


constituents, both public and private. 
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APPENDIX A 


APL programs used in chapter ITI. 


a) PROMATB 

The function PROMATB performs the calculations of 
ExanVlC——tidtmlco wm tetakeS aS input the billet structure 
B, the rotation rates R and the vector of withdrawal rates 
w. The unique promotion matrix Q which satisfies the basic 
balance equations and assumption A2;4 is then calculated. 


The function uses the function QCALC as a subroutine. 


V PROMATB5M;RO;ITQINV;G;K53V;RQ 


[1] "ARE THERE ANY CHANGES TO BILLET STRUCTURE?' 
[2] (OEM Meets: 
[3] L2:'SPECIFY LOCATION, PAY GRADE, NEW BILLET SIZE’ 
[4] FPG VIZ)ieV is) ,0e.Vel 
[5] "ANY MORE CHANGES?! 
[6] Reel) ) PLZ 
7 | Sa B xR 
[8] CLGINIEC 
i] HOUNVS @B NC (G,G) pl, Ged) -Q) [1s] 
[10] ‘CONTINUE?! 
[11] mae = |] J pLZ 
V 


Thee tfegmeorechances| 1m the billetsstructure on the 
Pi OlOmlOnm Grice mrclchy i Gigi nstanacemtGmecxanples 2aand 3. 
Came be ebsicrved sthrougn applreatiion ei PROMATB.) Line 1 
askse¥the user whether or not there are to be any changes in 
the billet Structure. As with other programs described in 
this section, an answer of 1 means affirmative and an answer 
ot Ul} means! a nerpative reply, line 3 tells the user who has 


indicated that changes are to be made to specify the changes. 


da 


-” 





With this information, a new Q is calculated, assuming that 


all other factors remain fixed. 


2) RELOCATE 
The function RELOCATE calculates the new billet structure 
pation ¢ a change in a particular element of B under the 
conditions of example 4. Recall that in Ss catieaile A@ot tchap- 
ter II, the billets in each grade were fixed, and thus the 
total population of the institution was fixed. The change 
in billets at one location for a particular grade was to be 
absorbed proportionately throughout the rest of the system. 
Following application of RELOCATE, the new billet struc- 


ture thus obtained is used as an input to PROMATB in order 


to generate a new promotion matrix Q. 


VeeEEOcATE: Gti: T:A 


[1] _ ‘WHAT LOCATION, PAY GRADE AND NEW BILLET SIZE?' 
[2] G+T[2],00, Ler en 00,T+[] 
[SI B[;G]<B{ ;G]x(A-T[3])+(A«(+/B) [G]) -B[L;G] 
[4] BiL;GiKT[3] 
V 


3) ROTCOMP 

Ties tuncesemmcn COMP calculateswehe new rotation rates 
[ORM anelCUlaRmearadcawhen themmotation rabettor personnel 
Of that grade at a specified location 1s changed. The rates 
Ger personneleet that gerade) at other location shift propor- 
Maciel lneron1sS No cC£feCt en the other parameters of 


the system. This was used in the calculations of example 8. 
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V ROLCOMP YG ib GA;C {BETA 


[1] ‘SPECIFY LOCATION GRADE AND NEW ROTATION RATE! 
[2] Cae 0 emo alc 0p, v (| 
[3] Aa RL GJEVi3)]-RML: Gl 
[4] BETA*+-AxC+(+/BxR) [G] -C<B[L;G]xR[L;G] 
[5] R[;G]«<R[;G]x1+BETA 
[6] R[L;G]<V[3] 
| V 
4)  PROMATR 


As opposed to PROMATB, which takes the billet structure 
an an input, PROMATR Gaxesm thle sLequluTLemnentsS mwerssevetal 
periods and calculates the promotion schemes for each of the 
periods. It is implicit in PROMATR that the tour lengths 
specified and the requirements are compatible, i.e. 
r(ai;k3;n) = r(i;k3;n-T(i;k)) for all locations i, grades k 
and time periods n. Thus PROMATR does the calculations 


under assumptions of weak stationarity. 


VPROMATR;V;IN;G;N;RQ;IQINV 

IN+1 

N<V[1] ,0p,G+V[3],0p,V+pRQN 

"ARE THERE ANY CHANGES TO SPECIFIED REQUIREMENTS? ' 
(OF) )eL2Z 
ibe UOPECIEY §BERTOD,LOCATION, PAY GRADE AND NEW REQUIRENEY 
RON TV LE] 5V1215V [3] J*vV (4) 00. V1 

PANY*®MORE CHANGES?! 

ra] Jolt 
L2:RQ+RON[IN; ; ] 

QCALC 

EXN[IN;;]+-RQ-A+RQ+.xQ 

ONT IN; ; ]<Q 

Reon (UNS lity ROre xW 

+(N>IN<IN+1) pL2 


MPEP rO ONAN BWNH 
IAN DRO RO ee Ss ss Fs ES ee Ss 


Re a. a nC 8 ee eet 
— SS | 


5) STATPRO 
STATPRO develops a stationary promotion scheme for systems 
TPM Mivemonlysweak Stationarity. The billet sSittructure sup- 


ported by a rotation scheme and requiyements as used in PROMATR 
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is used in PROMATB to generate the promotion scheme. The 
disparities between availabilities and requirements arising 
due to the forcing of a stationary promotion scheme onto a 


non-stationary structure are provided as an output. 


VeeoTATRBROS IN: EN:N:V 


[1] PROMATB 

[2] IN«1,09 ,N«V[1],00,V«oRQN 
[3] L1:RQ+RQN[IN; ;] 

[ 4 | EON RIN; ; |<-RO-AEROT RW 
[Sy +(N>IN+IN+1)oL1 


V 


6) QCALC 
QCALC is used as a subroutine in PROMATB and PROMATR, 


performing the actual calculations leading to the determina- 


tion of the promotion matrix Q. 


WEROCALC;K ;ROsG 
] G+ op IQINV<(+/RQ) = RO«++/RQ+. XW 
| emma. G) 00 ,0p,KezZ 
| aie ae OS a Oo Oa | er Na ie 
] Peer S K+) |e (OUR SKIN (KG 0o OTK Ki ai -1QINY 
[K-1]xQ[K-1;K]=IQINV[K] 
5] >+((G-1)>K+K+1)oL1 
6 | Q[G;G]+«1-W[G] 
V 
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